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ABSTRACT 

We present statistics of 133 faint 1.2-mm continuum sources detected in about 120 deep ALMA 
pointing data that include all the archival deep data available by 2015 June. We derive number 
counts of 1.2 mm continuum sources down to 0.02 mJy partly with the assistance of gravitational 
lensing, and find that the total integrated 1.2 mm flux of the securely identified sources is 22.9/g’g Jy 
deg -2 that corresponds to 104^25% of the extragalactic background light (EBL) measured by COBE 
observations. These results suggest that the major 1.2 mm EBL contributors are sources with 0.02 
mJy, and that very faint 1.2 mm sources with <0.02 mJy contribute negligibly to the EBL with the 
possible flattening and/or truncation of number counts in this very faint flux regime. To understand 
the physical origin of our faint ALMA sources, we measure the galaxy bias b g by the counts-in-cells 
technique, and place a stringent upper limit of b g < 3.5 that is not similar to b g values of massive DRGs 
and SMGs but comparable to those of UV-bright sBzKs and LBGs. Moreover, in optical and near- 
infrared (NIR) deep fields, we identify optical-NIR counterparts for 59% of our faint ALMA sources, 
majority of which have luminosities, colors, and the IRX-/3 relation same as sBzKs and LBGs. We 
thus conclude that about a half of our faint ALMA sources are dust-poor high-z galaxies as known as 
sBzKs and LBGs in optical studies, and that these faint ALMA sources are not miniature (U)LIRGs 
simply scaled down with the infrared brightness. 
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1. INTRODUCTION 


Since the infrared (IR) extragalactic background light 
(EBL) was first identified by the Cosmic Background Ex¬ 
plorer (COBE) satellite, it has been known that the to¬ 
tal en ergy of IR EBL is comparable to that of the optical 
EBL ( Puget et al.lll996HFixsen et al.lll998l:lHauser et al.l 


119981 Hauser fc Dwekl 120011 [Dole et al.l 120061 ). Because 
dusty star-forming galaxies at high redshift should sig¬ 
nificantly contribute to the IR EBL, the IR EBL is key 
to understanding the dusty-side of the cosmic star for¬ 
mation history and to constraining p hysical parameters 
of galaxy formation models (e.g., Granato et al.1 12004 
Baugh et alJ 120051: iFontanot et al.l 120071 iShimizu et al l 


2012t lHavward et al.ll2013D~ 


To reveal the origin of the IR EBL, we need reli¬ 
able number counts from bright to faint flux limits. 
The IR EBL of the low-frequency end, at submillime- 
ter (subnnn) and millimeter (mm) wavelengths, is one 
of the most advantageous wavelength regimes, due to 
the negative k-correction. Previous blank field obser¬ 
vations with single dish telescopes have revealed that 
the brig ht dusty high-z objec ts of submillimeter galaxies 
(SMGs; lLagache et al.ll2005l) are not major contributors 
of the EBL. SCUBA and LABOCA observations have 
resol ved 20 — 40% of t he EBL at submm wavelength 
(e.g., lEales et all 120001 iSmail et~aRl 120021 iGonoin et, al.l 
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120061 iKnudsen et al.ll2008l : iWeifi et aLll2009fl . Similar re¬ 
sults have been obtained by AzTEC observations that 
they have re solved 10 — 20% of the EBL at mm wave¬ 
length fe.g iPerera et al.l 120081 lHatsukade et al.l 120111 
iScott et al.ll2012D . These observations suggest that there 
should exist high-z populations different from SMGs, and 
that such populations are major EBL contributors in 
submm and mm wavelengths. 

There are two observational challenges for fully re¬ 
solving the EBL. First is the spatial resolution. The 
poor spatial resolution of the sing le-dish o bserva tions 
causes the serious source confusion (lCondonlll974l) that 
skews the observed number counts. Second is the 
sensitivity. Most of the submm and mm telescopes 
reach the limiting fluxes of ~ 1 mJy that corre¬ 
spond to the flux range of the SMGs. To detect 


submm sourc es fainter than ~ 1 mJy, there are many 
efforts (e.g.. ISmail et al.1 120021 [Chapman et al.1 120021 


Knudsen et al.l 20081 

Johansson et al. 2011: 

Chen et al.1 

20131). IKnudsen et al. 

(I2008D and Chen et al. 

(201311 ob- 


serve massive galaxy clusters at 850 /im for gravita¬ 
tionally lensed sources, and resolved 50%—180% of the 
EBL down to the intrinsic flux of S^o > 0.1 mJy. 
These results indicate that there exist faint (< 1 mJy) 
populations signi f icantly contributing to the EBL (e.g., 
ISmail et~aI1l2002t IKnudsen et al .1 [20081) . However, large 
uncertainties of the measurements still remain, due to 
the source blending and cosmic variance. Moreover, the 
physical origin of these faint submm and mm sources is 
unknown. 

The Atacama Large Millimeter / submillimeter Array 
(ALMA) enables us to investigate the faint submm and 
mm sources with the negligibly small uncertainties of the 
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source confusion and blending, having its capabilities of 
high angular resolution and sensitivity. The first ALMA 
mm search for faint sources below Simm ~ 1.0 mJy with 
nojensing effect has been conducted bv lHatsukade et all 
(2013). The study reports that they resolve ~80% of 
the EBL down to Si mm = 0.1 mJy in 20 targets re¬ 
siding in one blank f ield of the Subaru/X MM-Newton 
Deep Survey (SXDS; iFurusawa et al.ll2008l) . Due to the 
single field observations, this result would include un¬ 
known effects of cosmic variance. Subsequently, using 
multi-field deep ALMA 1mm maps of Ban ds 6 and 7 
down to th e flux limit of 5i m m = 0.1 mJy, lOno et ahl 
(12014H and iCarniani et al.l (|2015f) have claimed an EBL 
resolved fraction of ~60 % that is smaller than the one of 
lHatsukade et al.l (j2013f ). The uncertainties from cosmic 
variance are probably reduced with the multi-field data. 
Including the effects of the cosmic variance, only about 
a half of the EBL has been resolved with the present 
ALMA flux limit of Si mm = 0.1 mJy. 

There is another remaining issue about the IR EBL. 
Although these studies newly identify faint ALMA 
sources (Si mm <1.0 mJy) contributing to the half of 
the IR EBL, these studies do not clearly answer to 
the question about the connection between the faint 
ALMA sources and optically selected high-,? galaxies. 
Recent observations with Herschel have revealed that 
typical UV-selected galaxies such as Lyman-break galax¬ 
ies (LBGs) have a median total (8-1000 n m ) luminos¬ 
ity of Ljr ~ 2.2 x 10 n L w (|Reddv et al.l 1201*3 : iLee et all 
l2012t iDavies et ahll2013D . The stacking analysis of Her¬ 
schel and ALMA data has also shown that AT-selected 
galaxies, including star-forming BzK galaxies (sBzKs), 
have total IR luminosities of Air = (5 — 11) x 10 11 A 0 
(iDecarli et all [2014 ). These ranges of the IR luminosi¬ 
ties of LBGs and sBzKs correspond to the mm flux of 
S'imm ~ 0.1 — 1 mJy if we assume a modified black body 
with /3d = 1.8, dust temperature Id = 35K, and a source 
redshift z = 2.5. This mm flux range is similar to the one 
of the faint ALMA sources, implying that these sources 
could be mm counterparts of the optically selected galax¬ 
ies. 

There is another approach for characterizing the faint 
ALMA sources. Because the spatial distribution of galax¬ 
ies is related to the underlying distribution of dark 
haloes in the standard scenario of structure formation 
in ACDM universe, the clustering analysis is a power¬ 
ful tool for understanding the connection between vari¬ 
ous galaxy populations. Clustering analyses have been 
carried out over the past decade for bright SMGs ob¬ 
served with single-dish telescopes, and concluded that 
the clustering amplitudes, or galaxy biases, are large, 
and that the hosting dark halo ma sses of the b right 


SMGs are estimated to be ~ 10 ia M^ dWebb et al. 

2003; 

Blain et al. 2004 iWeifi et al.l 2009t 

Williams et al. 

2011: 

Hickox et al.ll2012l cf. iMillcr et all 

2015T). On the other 


hand, however, the dark halo properties of the faint 
A LMA sources a re poorly known, except for the result of 
lOno et aht (12014 ) who have obtained, for the first time, a 
meaningful constraint on a galaxy bias, b g < 4, for faint 
ALMA sources. In this way, the faint ALMA sources 
are not well studied. To understand the physical origin 
of faint ALMA sources, one should study faint ALMA 
sources with a more complete data set on the basis of 


individual sources as well as statistics. 

In this paper, we analyze the large dataset of multi- 
field deep ALMA data including the complete deep 
datasets archived by 2015 June. Sixty-six independent- 
field and one cluster maps are taken by ~ 120 ALMA 
pointings in Band 6/7 to reveal the origins of the EBL 
and faint ALMA sources in mm wavelength. Note that 
this is the first step for statistically investigating faint 
sources behind a lensing cluster with ALMA. The struc¬ 
ture of this paper is as follows. In Section [2 we describe 
the observations and the data reduction. Section [3] out¬ 
lines the method of the source extraction, our simula¬ 
tions for deriving the number counts, and the mass model 
development for the cluster. We compare the number 
counts from our and previous studies, and estimate the 
resolved fraction of the EBL at 1.2 mm in section HI In 
Section O we report the results of the clustering analysis 
for our faint ALMA sources. We investigate optical-NIR 
counterparts of our faint ALMA sources in Section [6] 
The summary of this study is presented in Section [3 

Throughout this paper, we assume a flat universe with 
fl m = 0.3, Da = 0.7, a s = 0.8, and H 0 = 70 km s _1 
Mpc -1 . We use magnitudes in the AB system (Oke & 
Gunn 1983). 

2. DATA AND REDUCTION 

We use 67 continuum maps obtained by ~ 120 pointing 
of the ALMA cycle 0-2 observations in Band 6/7 that ac¬ 
complish high sensitivities and angular resolutions. Ta¬ 
bles Q] and [2] present the summary and the detailed prop¬ 
erties, respectively, for these 67 continuum maps. We de¬ 
fine sub-datasets by the mapping modes and the depths 
because these two conditions would make different sys- 
tematics. For the mapping mode definitions, there are 
maps targeting field regions by single pointing observa¬ 
tions, referred to as ’field’ data, 4 and 62 out of which are 
taken from our programs and the ALMA archive, respec¬ 
tively. One map for a galaxy cluster taken by mapping 
observations, referred to as ’cluster’ data, is from the 
ALMA archive. 

For the depth definitions, we divide these 66 maps of 
’field’ data into two sub-datasets with low (< 60 pJy) and 
high (> 60 /xJy) noise levels, which we refer to as ’deep’ 
and ’medium-deep’ data, respectively. This is because 
the short-integration data would have more chances to 
contain systematic noise. As shown in Table[2] we use the 
sub-dataset names of A and B for the deep and medium- 
deep field data, respectively. The sub-dataset name of 
the deep cluster data is C. 

In this section, we describe the details of these 67 
continuum maps of ’deep’, ’medium-deep’, and ’cluster’ 
data, and present our reduction procedures. 

2.1. Our Data 

The 4 maps of the field data were obtained by our 
ALMA single pointing observations. Two out of the four 
maps were taken in the ALMA cycle-0 Band 6 obser¬ 
vations for the spectroscopicall y confirmed Lyq emitter 
(LAE) at z = 6.595, Himiko (lOuchi e t al.ll201 M and an 
LAE at z = 6.511 NB921-N- 7 9144 (lOuchi et alfeoTol PI 
R. Momose). See lOno et all (12014 for the summary of 
these observations. We also utilize two maps of newly ob¬ 
tained (ALMA cycle 2) field data of Bands 6 and 7 taken 
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Table 1 

Data Summary 


Sub-Dataset 

Number of Maps 

(i) 

(2) 

Deep data (A) 

41 

Medium-deep data (B) 

25 

Cluster data (C) 

1 


Note. — The sub-dataset names of A, 

B, and C are written in parenthesis. (1): 

ALMA maps with low (< 60 /iJy) and high (> 

60 AtJy) noise levels are referred to as Deep and 
Medium-deep data, respectively. (2): Number 
of the ALMA maps in each sub-dataset. 

for spect roscopica lly c onfirm ed LAEs at z = 5.7 (NB816- 
S-612 69: IQnc.hi ot a! .1 l200l and at 2 = 7.3 (NB101-S- 
2904; iKonno et al.1 12014 Ouchi et al. in preparation). 
The NB816-S-61269 observations were carried out in the 
ALMA program of #2012.1.00602.S (PI: R. Momose) on 
2014 May 20, Jun 19, and July 7 with 43 12-m anten¬ 
nae array in the range of 18 — 650 m baseline. The 
full width at half maximum (FWHM) of the primary 
beam was 20'/6. The available 7.5 GHz bandwidth with 
four spectral windows was centered at an observed fre¬ 
quency of 290.8 GHz (i.e. ~ 1.03 mm). J0217+0144 
and J2258-279 were observed as a flux calibrator, while 
J0006-0623 was used for bandpass calibrator. Phase cal¬ 
ibration was generally performed by using observations 
of J0217+0144 and J0215-0222. The total on-source ob¬ 
served time was ~ 2.3 hours. The NB101-S-2904 data 
were taken in the ALMA program of #2012.1.00088.S 
(PI: M. Ouchi) on 2014 July 22, August 6, 7, 14, and 
18, with 55 12-m antennae array in the extended con¬ 
figuration of 18-1300 m baseline. The full width at half 
maximum (FWHM) of the primary beam was 25'/3. The 
center of the observed frequency is 237.9 GHz (i.e. ~ 1.26 
mm). J0238+166, J2258-279, and J0334-401 were used 
for primary flux calibrators. Bandpass and phase cal¬ 
ibrations were performed with J0241-0815/J0238+1636 
and J0215-0222/J0217+0144, respectively. The total on- 
source integration time was ~ 5.2 hours. 

2.2. Archival Data 

To increase the number of ALMA sources, we make full 
use of ALMA archival data of cycles 0 and 1 that became 
public by 2015 June. We select all of the available ALMA 
Band 6/7 data sets fulfilling the following three criteria: 
1) The Id noise level is < 100 /iJy/beam for the contin¬ 
uum maps, 2) No bright (> 50cr) sources like AGN, QSO 
and SMG are included in the maps, and 3) The galactic 
latitude is high enough to avoid Galactic objects. The 
second criterion is important to reduce the chance for se¬ 
lecting residual side lobes of the bright sources that might 
remain even after the CLEAN algorithm on Common As¬ 
tronomy Software Applications (CASA; iMcMullin et al.1 
l2007t) package. Moreover, the second criterion signifi¬ 
cantly reduces the systematic uncertainties from galaxy 
clustering because it is reported that there is a clear ex¬ 
cess of submm number counts around bri g ht sources such 
as AG Ns and SMGs (jSilva et al.1 120151 : iSimoson et, al.1 
l2015aT) 

2.2.1. Field Data 


On the basis of our complete archival data search 
up to 2015 June with these three criteria, we col¬ 
lect the 62 maps of the fiel d data taken by si ngle¬ 
point ing observations: 2 fro m iWillott e t al.1 (12 013H. 2 0 
from lHatsukad e et al.1 (1201311 . 1 from lOta et al.1 (I2014T) 

2 from lHatsuka.de et al.1 (120141). 1 from iBeraer et a. 1.1 
(120141) . 3 from ICarniani et al.1 (1201511 . 20 from ALMA 
#2012.1.00076, 2 from ALMA#2012.1.00676.S, 2 from 
ALMA #2012.1.00323, 3 from ALMA #2012.1.00536, 
4 from ALMA #2012.1.00934, and 2 from ALMA 
#2012.1.00610. The ALMA #2012.1.00076 data were 
taken by PI K. Scott for 20 IR galaxies at 2 = 0.25 — 0.65 
in January 2014 with Band 6. The total on-source inte¬ 
gration is ~ 20 — 40 min for each map. The ALMA 
#2012.1.00676.S observations were carried out by PI C. 
Willott for 2 quasars at z ~ 6 in the end of November 
2013 with Band 6. The total on-source integration is 
~ 40 min. The ALMA #2012.1.00323.S observations 
were performed by PI G. Popping for 2 star-forming 
galaxies in March 2014 with Band 6. The total on- 
source integration is ^7 min for each source. The ALMA 
#2012.1.00536.S data were taken by PI C. Martin for 3 
LAEs at z = 5.7 with Band 7. The total on-source in¬ 
tegration is ~ 20 — 40 min for each source. The ALMA 
#2012.1.00934.S observations were conducted by PI B. 
Phillip for 4 normal star-forming galaxies at z = 2.23 
in January 2014 with Band 6. The total on-source 
integration is ~ 6 min for each source. The ALMA 
#2012.1.00610.S observations were carried out by PI B. 
Posselt for 2 nearby neutron stars in April and Jun 2014 
with Band 6. The total on-source integrations are ~ 30 
and 60 mins for each source. 

2.2.2. Cluster Data 

One-cluster data were taken for Abell 1689 (A1689) 
in ALMA cycle 0 and 1 observations (PI: J. Richard) 
of ALMA #2011.0.00319.S and #2012.1.00261.S. These 
two programs of A1689 included two different frequency 
settings, high and low tuning modes. The central- 
observed frequencies were chosen at 231 and 222 GHz 
(i.e. ~ 1.30 and 1.35 mm) for the high and the low tun¬ 
ing modes, respectively. Thus, a total of four observing 
sets were conducted for A1689 through ALMA cycle 0 
and 1. Each observing set comprises 2 — 4 sub-blocks 
of observations in which the source integration time is 
~ 10 — 30 minutes to make a mosaic map covering the 
highly magnified area of A1689 by ~50 pointings. The 
first observation set was obtained on 2012 June 17-18, 
and July 3-4 with 20 antennas in the high-tuning mode. 

Mars and Titan were observed as flux calibrators. The 
bandpass and phase calibrations were performed using 
observations of 3C 279. The second observation set was 
taken on 2012 July 14,15, 28, August 9, and 2013 Jan¬ 
uary 1 with 24 antennas in the low-tuning mode. The 
available 7.5 GHz bandwidth with four spectral windows 
was centered at an observed frequency of 222 GHz (i.e. 
^1.35 mm). Titan was observed as a flux calibrator. 
Bandpass and phase calibrations were performed with 
3C 279. The third observation set was obtained on 2014 
March 11 and April 4 with 34 antennas in the low-tuning 
mode. Titan was observed as a flux calibrator. Bandpass 
and phase calibrations were carried out with J1256-0547. 
The fourth observing set was performed on 2014 March 
22 and April 4 with 37 antennas in the high-tuning mode. 
































4 


Fujimoto et al. 

Table 2 

ALMA Maps of this Study 


Map ID 

Target 

^obs 

J'obs (Band) 

& cont 

beam 

size 

^1.2mm/ 'S'obs 

Ref. 

Project ID 



(mm) 

(GHz) 

(^zjy beam -1 ) 

("x 

") 






(i) 

(2) 

(3) 

(4) 

(5) 

(6) 


Deep Data (A) 

i 

BDF-3299 

1.30 

227 (6) 

8.5 

0.79x0.58 

1.28 

(a),(b) 

2012. A.00040.S 










2012.1.00719.S 

2 

GRB090423 

1.35 

222 (6) 

10.3 

0.98x0.75 

1.45 

( c )>( d ) 

2012.1.00953.S 

3 

NB101-S-2904 

1.26 

238 (6) 

11.4 

0.43x0.33 

1.16 

(e) 

2012.1.00088.S 

4 

NB816-S-61269 

1.03 

291 (7) 

13.0 

0.45x0.42 

0.63 

(f) 

2012.1.00602.S 

5 

rxj0806 

1.31 

229 (6) 

16.6 

0.78x0.63 

1.31 

(g) 

2012.1.00610.S 

6 

Himiko 

1.16 

259 (6) 

17.0 

0.81x0.56 

0.90 

0O,(i) 

2011.0.00115.S 

7 

BDF-521 

1.30 

230 (6) 

17.7 

0.68x0.51 

1.28 

(a),(b) 

2012.1.00719.S 

8 

IOK-1 

1.29 

232 (6) 

18.8 

1.08x0.78 

1.25 

(j),(k) 

2011.0.00767.S 

9 

rxj2143 

1.26 

238 (6) 

20.8 

0.60x0.43 

1.16 

(1) 

2012.1.00610.S 

10 

CFHQSJ2329-0301 

1.20 

250 (6) 

21.0 

0.73x0.63 

1.00 

( m )>( n ) 

2011.0.00243.S 

11 

SDF-46975 

1.22 

244 (6) 

21.2 

1.23x0.98 

1.05 

(°):( b ) 

2012.1.00719.S 

12 

ID239 

1.29 

232 (6) 

21.7 

1.88x0.90 

1.25 


2012.1.00076.S 

13 

GRB020819B 

1.18 

254 (6) 

22.6 

0.91x0.72 

0.95 

(P) 

2011.0.00232.S 

14 

CFHQSJ0055+0146 

1.09 

273 (6) 

23.3 

0.69x0.55 

0.75 

(m) 

2012.1.00676.S 

15 

ID182 

1.23 

244 (6) 

23.8 

1.33x0.89 

1.08 


2012.1.00076.S 

16 

ID220 

1.23 

244 (6) 

24.2 

1.33x0.88 

1.08 


2012.1.00076.S 

17 

ID244 

1.29 

232 (6) 

24.1 

1.87x0.90 

1.25 


2012.1.00076.S 

18 

CFHQSJ2229+1457 

1.12 

266 (6) 

24.9 

0.79x0.70 

0.81 

(m) 

2012.1.00676.S 

19 

ID247 

1.29 

232 (6) 

25.1 

1.87x0.90 

1.25 


2012.1.00076.S 

20 

ID209 

1.29 

232 (6) 

25.8 

1.23x1.01 

1.25 


2012.1.00076.S 

21 

CFHQSJ0210-0456 

1.20 

249 (6) 

27.2 

0.92x0.61 

1.00 

( m ),(n) 

2011.0.00243.S 

22 

MSDM80+3 

1.06 

283 (7) 

28.4 

0.84x0.63 

0.69 

(q) 

2012.1.00536.S 

23 

GRB051022 

1.15 

262 (6) 

28.6 

1.15x0.77 

0.88 

(p) 

2011.0.00232.S 

24 

ID217 

1.23 

244 (6) 

30.3 

1.33x0.93 

1.08 


2012.1.00076.S 

25 

ID225 

1.23 

244 (6) 

31.6 

1.34x0.93 

1.08 


2012.1.00076.S 

26 

SXDS19723 

1.27 

236 (6) 

38.5 

0.83x0.68 

1.19 

(r),(s) 

2011.0.00648.S 

27 

SXDS28019 

1.27 

236 (6) 

38.9 

0.84x0.68 

1.19 

(r),(B) 

2011.0.00648.S 

28 

SXDS22198 

1.27 

236 (6) 

39.0 

0.86x0.68 

1.19 

(r),(B) 

2011.0.00648.S 

29 

SXDS35572 

1.27 

236 (6) 

39.3 

0.84x0.68 

1.19 

(r),(s) 

2011.0.00648.S 

30 

SXDS103139 

1.27 

236 (6) 

39.7 

0.83x0.68 

1.19 

(r),( s ) 

2011.0.00648.S 

31 

NB921-N-79144 

1.22 

245 (6) 

50.1 

0.76x0.63 

1.05 

00,0) 

2012.1.00602.S 

32 

SXDS42087 

1.30 

231 (6) 

52.9 

0.80x0.66 

1.28 

(r),(s) 

2011.0.00648.S 

33 

SXDS79307 

1.30 

231 (6) 

53.7 

0.80x0.66 

1.28 

(r),(B) 

2011.0.00648.S 

34 

SXDS31189 

1.30 

231 (6) 

54.4 

0.81x0.66 

1.28 

(r),(B) 

2011.0.00648.S 

35 

MSDM29.5-5 

1.06 

283 (7) 

54.6 

0.71x0.53 

0.69 

(q) 

2012.1.00536.S 

36 

HiZELS-UDS-NBK-8806 

1.39 

216 (6) 

55.0 

1.86x0.89 

1.58 


2012.1.00934.S 

37 

HiZELS-UDS-NBK-11473 

1.39 

216 (6) 

55.3 

1.86x0.89 

1.58 


2012.1.00934.S 

38 

MSDM71-5 

1.06 

283 (7) 

56.1 

0.72x0.53 

0.69 

(q) 

2012.1.00536.S 

39 

HiZELS-UDS-NBK-13486 

1.39 

216 (6) 

56.4 

1.86x0.89 

1.58 


2012.1.00934.S 

40 

HiZELS-UDS-NBK-11961 

1.39 

216 (6) 

58.3 

1.86x0.89 

1.58 


2012.1.00934.S 

41 

ID125 

1.09 

273 (6) 

59.1 

1.30x0.73 

0.75 


2012.1.00076.S 

Medium-Deep Data (B) 

42 

ID 098 

1.09 

273 (6) 

60.1 

1.30x 

0.72 

0.75 


2012.1.00076.S 

43 

ID093 

1.09 

273 (6) 

60.3 

1.30x 

0.72 

0.75 


2012.1.00076.S 

44 

ID143 

1.09 

273 (6) 

61.7 

1.31x 

0.72 

0.75 


2012.1.00076.S 

45 

SXDS33244 

1.25 

240 (6) 

67.5 

1.04x 

0.64 

1.13 

00,0) 

2011.0.00648.S 

46 

ID163 

1.19 

252 (6) 

68.2 

1.36x 

0.92 

0.97 


2012.1.00076.S 

47 

ID204 

1.19 

252 (6) 

69.2 

L24x 

0.86 

0.97 


2012.1.00076.S 

48 

SXDS13316 

1.23 

244 (6) 

69.3 

0.87x 

0.64 

1.08 

(r),(s) 

2011.0.00648.S 

49 

ID158 

1.19 

252 (6) 

69.5 

1.35x 

0.92 

0.97 


2012.1.00076.S 

50 

ID192 

1.19 

252 (6) 

69.7 

L24x 

0.86 

0.97 


2012.1.00076.S 

51 

SXDS59863 

1.23 

244 (6) 

71.6 

0.85x 

0.65 

1.08 

(r),(s) 

2011.0.00648.S 

52 

ID 177 

1.13 

264 (6) 

73.0 

1.24x 

0.86 

0.83 


2012.1.00076.S 

53 

ID107 

1.13 

264 (6) 

73.8 

L24x 

0.86 

0.83 


2012.1.00076.S 

54 

SXDS67002 

1.23 

244 (6) 

74.0 

0.88x 

0.64 

1.08 

(r),(s) 

2011.0.00648.S 

55 

ID112 

1.13 

264 (6) 

75.4 

1.24x 

0.86 

0.83 


2012.1.00076.S 

56 

ID117 

1.13 

264 (6) 

75.5 

1.24x 

0.86 

0.83 


2012.1.00076.S 

57 

SXDS9364 

1.23 

244 (6) 

76.7 

0.87x 

0.64 

1.08 

(r),(s) 

2011.0.00648.S 

58 

SXDS13015 

1.23 

244 (6) 

81.3 

0.85x 

0.65 

1.08 

(r),(s) 

2011.0.00648.S 

59 

113083 

1.32 

226 (6) 

82.1 

1.18x 

0.72 

1.36 


2012.1.00323.S 

60 

SXDS68849 

1.25 

240 (6) 

82.3 

0.96x 

0.65 

1.13 

(r),(B) 

2011.0.00648.S 

61 

SXDS79518 

1.25 

240 (6) 

83.0 

1.02x 

0.65 

1.13 

(r),(s) 

2011.0.00648.S 

62 

1374240 

1.32 

226 (6) 

84.6 

1.20x 

0.72 

1.36 


2012.1.00323.S 

63 

SXDS101746 

1.25 

240 (6) 

86.7 

1.13X 

0.65 

1.13 

(r),(B) 

2011.0.00648.S 

64 

SXDS110465 

1.25 

240 (6) 

87.6 

1.13x 

0.64 

1.13 

(r),(s) 

2011.0.00648.S 

65 

SXDS1723 

1.25 

240 (6) 

94.0 

1.29x 

0.64 

1.13 

(r),(s) 

2011.0.00648.S 

66 

SXDS59914 

1.25 

240 (6) 

95.2 

1.20x 

0.65 

1.13 

0),(b) 

2011.0.00648.S 

Cluster Data (C) 

67 

A1689 

1.32 

231 (6) 

38.4-40.61 

0.96x0.74 

1.36 

(t) 

2011.0.00319.S 










2012.1.00261.S 

Note. — 

(1): Wavelength in the observed frame. (2): Frequency in the observed frame. (3): One sigma noise measured in each map 

before primary beam correction. (4): Synthesized beam size of our ALMA maps (weighting 

= ’natural’). (5): 

Ratio of the flux density at 

1.2 mm, Si. 

2 mm? to the one at the observed wavelength, S^bs? that is estimated with the modified blackbody spectrum of (3^ = 1.8, a dust 

temperature ol Ta = 35K, and a redshilt ot z = 2.5. (6): Reference, (a) Vanzella et al. 

2011); (b)IMaiolino et 

all (20151); (c)ITanvir et al. 

(2012^ (d) 

Berger et al. (20141); (e) IKonno et a 

l. (20141); (f) Ouchi et al. 420101); (g 

Haberl et al. (19981); (h) Ouchi 

et all (20131): (i) 

Ono et al. ( 

20141): (4) live et aU (2006): 

(k) Ota 

et al. (2014); 

(1) Zampieri et 

all (2001): 

(m) Willott et alJ ( 2010|); 

(n) Willott et al. 


(2013); (o) Ouchi ct al. (2009); (p) Hatsukade et al. (2014); (q) Martin ct al. (2006); (r) Yabc ct al. (2012); (s) Hatsukadc ct al. (2013); 
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Ceres was observed as a flux calibrator. Bandpass and 
phase calibrations were done with J1256-0547. Through¬ 
out these four observing sets, the total on-source integra¬ 
tion is 5.5 hour. 

2.3. Data Reduction 

Basically, the data are reduced with CASA version 
4.3.0 in a standard manner with the scripts used for the 
data reduction provided by the ALMA observatory. In 
this process, we also use previous CASA versions from 
3.4 to 4.2.2, if we find problems in the final image that 
the noise level is significantly higher than the calibrated 
products provided by the ALMA observatory, or that 
there remain striped patterns. Exceptionally for the cy¬ 
cle 0 data with the high noise level or the patterned noise, 
we use the calibrated data produced by the ALMA obser¬ 
vatory. Similarly for the data of ALMA $=2011.0.00648.S 
(PI: Ohta), we use a re-calibrated data provided by Seko 
et al. because they find that the coordinate of a phase 
calibrator of #2011.0.00648.S is wrong (~ 0'/3), which 
causes positional offs ets (Seko et al. in preparation; 
lHatsukade et al.ll201 5<l. 

Our CASA reduction has three major steps: bad data 
flagging, bandpass calibration, and gain calibration with 
flux scaling. In the first step, we remove the data of 
the shadowed antennas and the edge channels of spec¬ 
tral windows. We also do not use the unreliable data 
such with jumps or a low phase/amplitude gain. We ap¬ 
ply the flaggings that are shown in the scripts, but no 
additional flaggings. This is because we find that the 
flaggings in the scripts are good enough for our scientific 
goals. The second step is the bandpass calibration. Af¬ 
ter we calibrate the phase time variation on the bandpass 
calibrator scan, we obtain the bandpass calibration in the 
phase and amplitude. The final step is the gain calibra¬ 
tion. We estimate the time variations of the phase and 
amplitude on the phase calibrators, and then transfer the 
calibration to the target source by applying the linear in¬ 
terpolation with the results of the phase calibrators. The 
water vapor radiometer is used for the correction of the 
short-time phase variations that cannot be traced by the 
phase calibrators. We estimate flux scaling factors, us¬ 
ing solar-system objects and bright QSOs with the flux 
models of ’Bulter-JPL-Horizons 2012’ and ALMA Cal¬ 
ibrator Source Catalog, respectively. The difference of 
the opacity between the flux calibrators and the target 
sources is corrected with the system noise temperature 
(T sys ) measurements. The systematic flux uncertainty is 
typically 10% in bands 6 and 7 (ALMA proposer’s guide 
0 ). 

We perform Fourier transformation for the original uv- 
data to create ’’dirty” maps. The continuum maps are 
made using all the line free channels of the four spectral 
windows. We process the maps with the auto CLEAN 
algorithm down to the depth of 3a noise levels of the dirty 
maps, using natural weighting. The final cleaned maps 
achieve angular resolutions of 0'/45 x 0'/33 —1"88 x 0'.'90. 
The sensitivities of these ALMA maps range from 8.5 to 
95.2 /ijy beam -1 before the primary beam corrections. 

Because the observations of the cluster data were taken 
place in 4 different epochs, we recalculate the data 

4 Section A7: https://almascience.nao.ac.jp/documents-and- 
tools / cycle -1 / alma-proposers-guide 



Figure 1. False-color image of A1689 cluster (red: 2775 , green: 
Vq 25 1 red: -B 475 ). The critical lines for background sources at 
2 : = 2.5 are shown with the magenta lines. The white curves are 
the contours indicating the 50% sensitivity of the deepest part in 
the mosaic map. The white thick lines and ticks indicate the 9 
regions of our noise estimates. 

weights with the statwt task in CASA based on its vis¬ 
ibility scatters which include the effects of integration 
time, channel width, and systematic temperature. After 
applying the data weights to the datasets, we combine 
these two datasets with the concat task in CASA. Since 
the mosaic mapping observations were conducted for the 
cluster, the noise levels vary by positions. To estimate 
the noise levels of the cluster data, we divide the cluster’s 
mosaic map in 9 regions as shown in Figure [I] The data 
of the 9 regions achieve the la noise levels in the range 
of 38.4-40.6 /My. 

3. DATA ANALYSIS 

We show our data analyses in the following subsec¬ 
tions. The wavelengths of our mm maps fall in the range 
of 1.03 — 1.32 mm. Our data set is composed of 63 Band 
6 maps and 4 Band 7 maps (Table 2), and the majority 
of our data are taken with Band 6. Because our ma¬ 
jor data of the ALMA Band 6 maps have the average 
wavelength of 1.23 mm, we derive the number counts 
with a flux density at 1.2 mm, 5i.2mm- We scale the 
flux densities of our mm maps to Si. 2 mm with the flux 
density ratios summarized in Tabic [2] For the data of 
the previous studies defined with the wavelengths dif¬ 
ferent from 1.2 mm, we estimate Si. 2 mm with the ra¬ 
tios of Si.2mm/Sg50itm = 0.37, Si.2mm/Ss70iim = 0.39, 
Sl.2mm/Sl.lmm — 0.77, and Si.2mm/Sl.3mm — 1.28. All 
of these flux density ratios are calculated based on a mod¬ 
ified blackbody whose spectral index /3d and dust temper- 
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ature Th are similar to th ose of typical SMGs; fig = 1.8 
(e.g., iChapin et al.ll2009l : Planck Collaboration 1201 111 
and T d = 35 K (e.g.. iKovkcs et all 120061 : iCoppin et alJ 
1200811 . Here we assume a source redshif t of z = 2.5 that 
is a median redshift value of SMGs ('e.g.-IChapman et al.l 
120051 lYun et aLll2012t ISimpson et al.ll2014fl . 

3.1. Source Detection 

We conduct source extraction for our ALMA maps 
before pri mary beam corrections with SExtractor ver¬ 
sion 2.5.0 (IBertin &; Arnoutsl 11 99611 . textcolorblue The 
source extraction is carried out in high sensitivity re¬ 
gions. For the single-pointing maps presented in Sections 
12.11 and 12.2.11 we use the regions with the primary beam 
sensitivity greater than 50%. For the mosaic map shown 
in Section 12.2.21 we perform source extraction where the 
relative sensitivity to the deepest part of the mosaic map 
is greater than 50%. 

We identify sources with a positive peak count above 
the 3.0cr level, assuming the sources are not resolved out 
in our ALMA maps. The catalog of these sources is re¬ 
ferred to as the 3.0a-detection catalog. From our 3.0a- 
detection catalog, we remove the objects located at the 
map centers that are main science targets of the archived- 
data ALMA observations. 


3.2. Spurious Sources and Source Catalog 

The 3.0a-detection catalog should include many spuri¬ 
ous sources, because the detection is simply determined 
with a peak pixel count. To evaluate the number of 
the spurious sources, we conduct negative peak analysis 
(lHatsukade et al.l 120131 : IQno et all 12014 ICarniani et alJ 
I2015H for the deep (A), medium-deep (B), and cluster 
(C) datasets. We multiply —1 to each pixel value of 
our ALMA maps, and perform the source extraction for 
negative peaks in the same manner as those in Section 
13.11 The numbers of negative and positive peaks in our 
ALMA maps are shown in Figure [5] The excess of the 
positive to negative peak numbers is regarded as the 
real source numbers. We model the negative and pos¬ 
itive peak distributions of Figure [2] with the functions of 
signal-to-noise ratio SNR, 

AT p = a n x i 0 - b " xSNR , (1) 

./Vpp = a n x 10- b ” xSNR + a p x io~ b P xSNR , (2) 

where N^ p and N pp are the cumulative numbers of nega¬ 
tive and positive peaks, respectively, and a n , b n , a p , and 
b p are free parameters. 

We estimate spurious source rates / sp defined by the 
ratio of negative to positive peak numbers, 


/sp(SNR) 


-^np(SNR) 

JV p ySNR)' 


(3) 


We evaluate / sp with the best-fit functions of eqs. m 
and ©, and present / sp in Figure [3] 

Figures [3] indicates that the field data of deep (A) and 
medium-deep (B) have an almost identical distribution 
of the spurious source rates, while the cluster data (C) 
shows spurious source rates higher than (A) and (B) at 
a given SNR. This result would suggest that spurious 
source rates do not depend on the data depth, but the 





Figure 2. Total numbers of positive and negative peaks in the 
deep (A), medium-deep (B), and cluster (C) sub-datasets as a func¬ 
tion of SNR. The thick solid and dashed lines show the total num¬ 
bers of the positive and negative peaks, respectively. The thin 
solid and dashed curves denote the best-fit model functions. The 
red lines are selection limits for our source catalog. See text for 
more details. 
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Figure 3. Cumulative spurious source rates as a function of SNR. 
The solid, dashed, and dotted lines indicate the deep (A), medium- 
deep (B), and cluster (C) sub-datasets, respectively. The red lines 
at SNRs of 3.4 and 3.9 are the selection limits (SL) of our source 
catalog. 

mapping modes. The complex distribution of the depths 
in the mapping data of (C) would have the relation of 
SNR and spurious source rate that is different from the 
smooth distribution of the depth in the single pointing 
data of (A) and (B). We use sources down to an SNR 
whose spurious source rate is < 70%, and adopt 3.4, 
3.4, and 3.9er levels for our selection limits of the data 
(A), (B), and (C), respectively. Applying these selection 
limits to our 3.0er-detection catalog, we obtain a source 
catalog consisting of 133 sources: 122 from the field data 
(A and B) and 11 from the cluster data (C). The source 
catalog of our 133 faint ALMA sources is presented in 
Table El 

Note that A1689-z Dl, an LBG at z = 7.5 behind 
A1689 (IWatson et al.ll2015l ). is not included in our source 
catalog because A1689-zDl is located outside of the pri¬ 
mary beam with a > 50% sensitivity that is one of our 
selection criteria. Nevertheless, we have looked at the 
position of Al689-zDl in our data, and found that there 
is a source with a n SNR of 4 w h ose flu x is the same as 
the one derived bv lWatson ct al.1 (12015 ) within the ~ ler 
uncertainty. We confirm the ALMA source detection at 
the position of Al689-zDl. 

3.3. Completeness and Flux Boosting 

We estimate the detection completeness by Monte- 
Carlo simulations. We put a flux-scaled synthesized 
beam into a map as an artificial source on a random posi¬ 
tion. These artificial sources have SNRs ranging from 3.0 
to 7.0 with an SNR step of 0.2 dex. If an artificial source 
is extracted within a distance of the synthesized beam 
size from the input position, we regard that the source 
is recovered. These completeness estimations are con¬ 
ducted in the data before primary beam corrections, be¬ 
cause we perform the source extractions in the real data 
(Section l3.il) uncorrected for primary beam attenuations. 
For all of our ALMA maps, we iterate this process 100 
times at an SNR bin. Figure [4] displays our completeness 



SNR 

Figure 4. Completeness as functions of SNR for three sub¬ 
datasets. The circles (red line), squares (blue line), and diamonds 
(green line) represent the completeness values (the best-fit func¬ 
tions) of the ID 6, 61 and 67 maps, respectively. 

estimates for three sub-datasets. We fit the function of 
1 — exp(aSNR — b) to our completeness estimates, where 
a and b are free parameters. The best-fit functions are 
also shown in Figure [4] 

We examine the flux boosting that is caused 
by the confusions of t he un detected faint sources 
(jAustermann ct ahl 12009 . I2010D . We use the detected 
objects of the Monte-Carlo simulations conducted for 
the completeness estimates, and measure the output flux 
densities at the input positions. We find that the ratios 
of the input to output flux densities are almost unity 
with a variance of only 5% over the wide SNR range of 
our ALMA maps. This small flux boosting is probably 
found, because the angular resolution of ALMA is high. 
We thus conclude that the effect of flux boosting is neg¬ 
ligibly small. 

3.4. Flux Measurement 

We estimate the source fluxes of our 133 objects 
from the integrated flux densities measured with the 
2D Gaussian-fitting routine of imfit in CASA. Because 
the integrated flux density values for low-SNR sources 
would include large systematic uncertainties originated 
from the spatially-correlated positive noise at around the 
source position, we use the integrated flux densities of 
imfit only for the sources which are identified as ’re¬ 
solved’ with an SNR of > 5. For the rest of the sources, 
we use the peak fluxes of the imfit best-fit Gaussian 
measurements. We confirm that the peak flux values 
of imfit , S'imfit, agree with those of SExtractor, Ssex, 
within the la error, S'imfit/S'sEx = 0.93 ± 0.15. 

Because the sources in the cluster data of A1689 would 
be distorted by the gravitational lensing effects, we make 
a low-resolution map of A1689 produced with a 1" cir¬ 
cular gaus sian uv-taper t hat is the same as the one pro¬ 
duced by IWatson et al.l (I2015D . The taped beam size 
is 1"31 x 1"12. In this map, we perform the imfit 
flux measurements (as described above) at the source 
positions determined in the original map. Here we de- 
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fine flux measurements in the tapered and non-tapered 
maps as .Staler and S/lo taper, respectively. We then ob¬ 
tain »Staper / S no taper = 1-01 ± 0.19 on average, if we omit 
two sources showing signatures of source confusions by 
neighboring objects. However, the strongly-lensed source 
of Al shows a systematically high fraction of Staper / 
iS'notaper = 1.23. Thus, we decide to use .Staper f° r the 
cluster data of A1689, except for those two confusion 
sources. We also make low-resolution maps for some of 
the field data, and perform the same test of tapering. Be¬ 
cause we find that .Staper and »S no taper are almost identical 
in these data, we use <S no taper for the field data. 

We test the reliability of our flux measurements, per¬ 
forming flux recovery simulations. First, we select the 
typical data sets of ID22 and ID55 that have beam sizes 
of ~ 0"8 and ~ 1"2 corresponding to the two peaks of bi- 
modal size distribution for our data. We then create 200 
input sources of elliptical Gaussian model profiles with 
the uniform distribution of SNRs in 3 — 7 and major- 
axis sizes of 0"2 and 0"5. Here, the input source sizes 
are chosen from rec ent A LMA high-reso l ution s tudies. 
Ilkarashi et al.l (|2014H and ISimpson et al.1 (|2015bH claim 
that the median sizes o f individual submm so urces are 
~ 0"2 — 0"3 in FWHM. ISimpson et al.l (|2015bH also find 


that the stacked image of submm sources with a low SNR 
of 4 — 5 has the size of FWHM 0"35io;i[o- Thus, we 
investigate the model sources of 0"2 and 0"5 that are 
taken from the median individual source size and the ltr 
upper limit stacked source size. We add these model 
sources to ID22 and ID55 maps at random positions. Fi¬ 
nally, we obtain flux densities of the model sources in 
the same manner as our flux measurements for the real 
sources, and compare these flux measurements with the 
input source fluxes. In the case of the 0"2 input source 
size, we find that the flux recovery ratio (defined by the 
ratio of the output to input fluxes) is 1.0 ±0.1 (1.0 ±0.2) 
for ID22 (ID55) data. Similarly, in the case of the 0"5 
input source size, we obtain 0.8 ±0.1 (0.9 ±0.2) for ID22 
(ID55) data. These results suggest that the flux recovery 
ratios are nearly unity within the 1 — 2er uncertainties, 
and that our flux measurements are reliable. 

Note that missing fluxes of interferometric observa¬ 
tions are negligible. ALMA proposer’s guide 0 shows 
that there exist missing flux effects for sources with 
~ 2"0 — 3"0 at Band 6/7 for the shortest-baseline con¬ 
figuration that we use. The source sizes of ~ 2'/0 — 3'/0 
are significantly larger than our ALMA sources whose 
angular sizes are < 0"5, and the missing fluxes should be 
negligible. 


Table 3 

Source Catalog 


ID 

Map ID 

R.A.(J2000) 

Decl.(J2000) 

gobs 

ocorr 
° 1.2 mm 

SNR 

/sp 

( 1 ) 

( 2 ) 



(mjy) 

(3) 

(mjy) 

(4) 

(5) 

( 6 ) 

1 

1 

337.049805 

-35.169643 

0.06 

± 

0.02 

0.07 

± 

0.02 

3.6 

0.57 

2 

1 

337.050232 

-35.164864 

0.41 

± 

0.01 

0.68 

± 

0.07 

39.3 

0.00 

3 

1 

337.050293 

-35.170044 

0.07 

± 

0.02 

0.09 

± 

0.02 

4.0 

0.34 

4 

1 

337.053314 

-35.165119 

0.06 

± 

0.01 

0.12 

± 

0.02 

5.7 

0.00 

5 

1 

337.054413 

-35.167576 

0.05 

± 

0.01 

0.07 

± 

0.02 

4.1 

0.27 

6 

1 

337.054718 

-35.166332 

0.05 

± 

0.01 

0.06 

± 

0.02 

3.5 

0.66 

7 

2 

148.887314 

18.150942 

0.05 

± 

0.01 

0.04 

± 

0.02 

3.6 

0.63 

8 

2 

148.891403 

18.151211 

0.06 

± 

0.02 

0.06 

± 

0.02 

3.4 

0.71 

9 

3 

34.446007 

-5.008159 

0.05 

± 

0.01 

0.05 

± 

0.02 

3.5 

0.70 

10 

3 

34.446648 

-5.007715 

0.05 

± 

0.01 

0.04 

± 

0.02 

3.4 

0.70 

lit 

4 

34.437332 

-5.492185 

0.06 

± 

0.02 

0.04 

± 

0.01 

3.7 

0.53 

12 + 

4 

34.437649 

-5.494778 

0.06 

± 

0.02 

0.03 

± 

0.01 

3.5 

0.67 

13+ 

4 

34.437805 

-5.494339 

0.06 

± 

0.01 

0.03 

± 

0.01 

3.8 

0.48 

14+ 

4 

34.439056 

-5.494986 

0.06 

± 

0.02 

0.03 

± 

0.01 

3.6 

0.60 

15+ 

4 

34.440075 

-5.494945 

0.07 

± 

0.02 

0.04 

± 

0.01 

3.6 

0.63 

16 

5 

121.593948 

-41.376675 

0.10 

± 

0.03 

0.08 

± 

0.04 

3.6 

0.60 

17 

5 

121.596870 

-41.375671 

0.06 

± 

0.02 

0.08 

± 

0.02 

3.6 

0.63 

18 

5 

121.599480 

-41.373417 

0.08 

± 

0.02 

0.10 

± 

0.03 

3.7 

0.55 

19 

5 

121.599564 

-41.375568 

0.07 

± 

0.02 

0.08 

± 

0.03 

3.6 

0.61 

20 

5 

121.600471 

-41.375629 

0.08 

± 

0.02 

0.09 

± 

0.03 

3.6 

0.60 

21 

7 

336.940674 

-35.119549 

0.11 

± 

0.03 

0.13 

± 

0.04 

3.8 

0.48 

22 

8 

201.001526 

27.413574 

0.10 

± 

0.03 

0.11 

± 

0.04 

3.5 

0.70 

23 

8 

201.001526 

27.413574 

0.10 

± 

0.03 

0.11 

± 

0.04 

3.5 

0.70 

24 

8 

201.001526 

27.413574 

0.10 

± 

0.03 

0.11 

± 

0.04 

3.5 

0.70 

25 

9 

325.760986 

6.904217 

0.13 

± 

0.04 

0.17 

± 

0.05 

3.5 

0.68 

26 

9 

325.762177 

6.902022 

0.14 

± 

0.04 

0.17 

± 

0.05 

3.6 

0.62 

27 

9 

325.763336 

6.905892 

0.08 

± 

0.02 

0.10 

± 

0.03 

3.4 

0.72 

28 

9 

325.767029 

6.902850 

0.15 

± 

0.04 

0.16 

± 

0.05 

3.6 

0.60 

29 

10 

352.281616 

-3.033988 

0.13 

± 

0.04 

0.12 

± 

0.04 

3.4 

0.71 

30 

10 

352.284851 

-3.031080 

0.16 

± 

0.03 

0.16 

± 

0.03 

5.9 

0.00 

31 

10 

352.285248 

-3.030138 

0.17 

± 

0.04 

0.17 

± 

0.04 

4.6 

0.08 

32 

10 

352.285309 

-3.035150 

0.11 

± 

0.03 

0.10 

± 

0.03 

3.8 

0.47 

33 

11 

200.929321 

27.340736 

0.15 

± 

0.04 

0.13 

± 

0.04 

3.9 

0.40 

34 

12 

150.547440 

1.910628 

0.13 

± 

0.04 

0.13 

± 

0.05 

3.7 

0.55 

35 

13 

351.831116 

6.267834 

0.13 

± 

0.03 

0.12 

± 

0.03 

4.2 

0.22 

36 

13 

351.833496 

6.267574 

0.16 

± 

0.04 

0.15 

± 

0.04 

3.7 

0.53 


5 Table A2 of the ALMA proposers’ guide: 
https: / / almascience. nao. ac. j p / documents-and-tools / cycle -1 / alma- 


proposers-guide 
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Table 3 — Continued 


ID 

Map ID 

R.A.(J2000) 

Decl.(J2000) 

gobs 

•^obs 

(mjy) 

ocorr 
° 1.2mm 

SNR 

fsp 





(mjy) 



(1) 

(2) 




(3) 



(4) 


(5) 

(6) 

37 

14 

13.760888 

1.771644 

0.10 

± 

0.03 

0.08 

± 

0.02 

4.0 

0.35 

38 

14 

13.761566 

1.768838 

0.16 

± 

0.04 

0.12 

± 

0.04 

3.7 

0.52 

39 

14 

13.763390 

1.770106 

0.11 

± 

0.03 

0.08 

± 

0.03 

3.7 

0.52 

40 

14 

13.764014 

1.770272 

0.13 

± 

0.03 

0.09 

± 

0.03 

3.7 

0.50 

41 

15 

149.775772 

1.736944 

0.16 

± 

0.03 

0.18 

± 

0.04 

5.1 

0.02 

42 

15 

149.775955 

1.738510 

0.15 

± 

0.04 

0.17 

± 

0.05 

3.9 

0.41 

43 

15 

149.776794 

1.736661 

0.09 

± 

0.03 

0.07 

± 

0.03 

3.6 

0.63 

44 

18 

337.254547 

14.950516 

0.17 

± 

0.05 

0.14 

± 

0.04 

3.8 

0.48 

45 

18 

337.256073 

14.954236 

0.19 

± 

0.03 

0.21 

± 

0.03 

6.1 

0.00 

46 

19 

149.722565 

2.606706 

0.11 

± 

0.03 

0.10 

± 

0.04 

3.4 

0.70 

47 

19 

149.724045 

2.602218 

0.21 

± 

0.05 

0.28 

± 

0.07 

4.1 

0.26 

48 

20 

149.792572 

2.103762 

0.12 

± 

0.03 

0.10 

± 

0.04 

4.0 

0.30 

49 

21 

32.553406 

-4.939113 

0.12 

± 

0.03 

0.11 

± 

0.03 

3.9 

0.37 

50 

21 

32.553787 

-4.939454 

0.11 

± 

0.03 

0.15 

± 

0.03 

3.7 

0.56 

51 

21 

32.555397 

-4.937203 

0.13 

± 

0.03 

0.13 

± 

0.04 

3.7 

0.54 

52t 

22 

150.127991 

2.288839 

0.13 

± 

0.04 

0.08 

± 

0.03 

3.4 

0.72 

53 

23 

359.017120 

19.607481 

0.10 

± 

0.03 

0.08 

± 

0.03 

3.4 

0.71 

54 

24 

150.406784 

2.034714 

0.14 

± 

0.04 

0.12 

± 

0.04 

3.5 

0.64 

55 

25 

150.393936 

1.998264 

0.15 

± 

0.04 

0.16 

± 

0.05 

3.7 

0.56 

56 

26 

34.100346 

-5.154171 

0.15 

± 

0.04 

0.14 

± 

0.05 

3.5 

0.64 

57 

26 

34.100658 

-5.156616 

0.16 

± 

0.05 

0.17 

± 

0.06 

3.5 

0.66 

58 

26 

34.102978 

-5.155593 

0.15 

± 

0.04 

0.15 

± 

0.05 

3.4 

0.70 

59 

27 

34.034069 

-5.104988 

0.17 

± 

0.04 

0.18 

± 

0.06 

3.8 

0.45 

60 

28 

34.472076 

-4.711596 

0.26 

± 

0.07 

0.27 

± 

0.09 

3.5 

0.65 

61 

28 

34.473747 

-4.713997 

0.16 

± 

0.04 

0.20 

± 

0.06 

3.6 

0.62 

62 

29 

34.392979 

-5.043759 

0.17 

± 

0.04 

0.19 

± 

0.06 

3.7 

0.51 

63 

29 

34.396526 

-5.042582 

0.21 

± 

0.06 

0.23 

± 

0.07 

3.6 

0.59 

64 

29 

34.396526 

-5.046781 

0.25 

± 

0.07 

0.26 

± 

0.09 

3.4 

0.72 

65 

31 

34.612003 

-4.588018 

0.25 

± 

0.07 

0.28 

± 

0.08 

3.5 

0.69 

66 

31 

34.612686 

-4.583063 

0.36 

± 

0.07 

0.43 

± 

0.09 

4.9 

0.04 

67 

31 

34.613739 

-4.588505 

0.31 

± 

0.09 

0.31 

± 

0.10 

3.5 

0.68 

68 

31 

34.613930 

-4.585264 

0.19 

± 

0.06 

0.18 

± 

0.06 

3.5 

0.70 

69 

33 

34.272991 

-4.859830 

0.29 

± 

0.08 

0.34 

± 

0.11 

3.5 

0.68 

70 

33 

34.274231 

-4.860201 

0.32 

± 

0.09 

0.43 

± 

0.12 

3.6 

0.59 

71 

33 

34.274761 

-4.859237 

0.26 

± 

0.07 

0.32 

± 

0.09 

3.8 

0.46 

72 

34 

34.305080 

-5.067344 

0.29 

± 

0.07 

0.39 

± 

0.10 

3.9 

0.41 

73 

34 

34.305538 

-5.069685 

0.23 

± 

0.06 

0.33 

± 

0.09 

3.6 

0.61 

74 

34 

34.305553 

-5.067293 

0.29 

± 

0.07 

0.22 

± 

0.09 

4.2 

0.19 

75 

34 

34.307556 

-5.071889 

0.33 

± 

0.09 

0.48 

± 

0.12 

3.6 

0.61 

76 

34 

34.309822 

-5.069085 

0.29 

± 

0.08 

0.23 

± 

0.11 

3.5 

0.65 

77+ 

35 

230.738800 

-0.125928 

0.35 

± 

0.10 

0.23 

± 

0.07 

3.5 

0.69 

00 

35 

230.739777 

-0.125576 

0.28 

± 

0.08 

0.20 

± 

0.06 

3.6 

0.59 

-t— 

Oi 

35 

230.743118 

-0.127994 

0.28 

± 

0.08 

0.18 

± 

0.06 

3.5 

0.66 

80 

36 

34.327038 

-5.131583 

0.21 

± 

0.06 

0.27 

± 

0.09 

3.8 

0.48 

81 

37 

34.401154 

-5.152167 

0.33 

± 

0.10 

0.37 

± 

0.15 

3.4 

0.72 

82 

37 

34.401531 

-5.154098 

0.30 

± 

0.08 

0.34 

± 

0.13 

3.7 

0.55 

83 

37 

34.403275 

-5.149917 

0.36 

± 

0.11 

0.61 

± 

0.18 

3.4 

0.71 

84+ 

38 

231.036255 

-0.176025 

0.39 

± 

0.11 

0.22 

± 

0.08 

3.5 

0.67 

85+ 

38 

231.037949 

-0.178142 

0.22 

± 

0.06 

0.15 

± 

0.04 

3.7 

0.54 

86 

40 

34.414261 

-5.201880 

0.34 

± 

0.10 

0.56 

± 

0.16 

3.5 

0.65 

87 

40 

34.419544 

-5.199385 

0.31 

± 

0.08 

0.52 

± 

0.14 

3.8 

0.46 

88 

41 

149.590530 

2.807018 

0.28 

± 

0.07 

0.23 

± 

0.05 

4.2 

0.21 

89 

41 

149.591949 

2.805352 

0.26 

± 

0.07 

0.21 

± 

0.05 

4.0 

0.31 

90 

42 

150.688171 

2.681118 

0.37 

± 

0.07 

0.26 

± 

0.06 

5.4 

0.01 

91 

45 

34.194702 

-5.056259 

0.43 

± 

0.12 

0.39 

± 

0.15 

3.5 

0.71 

92 

45 

34.197121 

-5.059778 

0.31 

± 

0.08 

0.37 

± 

0.10 

3.7 

0.55 

93 

46 

150.409714 

2.358208 

0.30 

± 

0.09 

0.32 

± 

0.09 

3.5 

0.71 

94 

47 

149.770050 

1.804996 

0.32 

± 

0.09 

0.31 

± 

0.10 

3.4 

0.72 

95 

48 

34.411190 

-4.745051 

0.30 

± 

0.08 

0.14 

± 

0.09 

3.8 

0.44 

96 

48 

34.411701 

-4.745893 

0.31 

± 

0.08 

0.34 

± 

0.09 

4.0 

0.31 

97 

49 

150.600494 

2.444365 

0.47 

± 

0.13 

0.53 

± 

0.14 

3.6 

0.59 

98 

51 

34.440269 

-4.912953 

0.39 

± 

0.11 

0.40 

± 

0.12 

3.6 

0.60 

99 

51 

34.441357 

-4.912045 

0.32 

± 

0.08 

0.29 

± 

0.09 

3.9 

0.43 

100 

51 

34.441635 

-4.913539 

0.46 

± 

0.13 

0.45 

± 

0.14 

3.6 

0.62 

101 

51 

34.442795 

-4.911066 

0.53 

± 

0.08 

0.56 

± 

0.11 

6.2 

0.00 

102 

51 

34.443920 

-4.909602 

0.41 

± 

0.12 

0.35 

± 

0.13 

3.5 

0.65 

103 

52 

150.151459 

1.751314 

0.49 

± 

0.14 

0.50 

± 

0.13 

3.6 

0.64 

104 

54 

34.759140 

-4.832398 

0.32 

± 

0.09 

0.37 

± 

0.11 

3.5 

0.71 

105 

55 

150.022736 

2.027032 

0.53 

± 

0.08 

0.46 

± 

0.08 

6.8 

0.00 

106 

57 

34.141144 

-5.232085 

0.47 

± 

0.13 

0.42 

± 

0.15 

3.6 

0.60 

107 

58 

34.303791 

-5.161064 

0.47 

± 

0.14 

0.58 

± 

0.16 

3.4 

0.73 

108 

59 

150.201645 

2.021740 

0.57 

± 

0.15 

0.72 

± 

0.22 

3.7 

0.53 

109 

60 

34.250214 

-4.804325 

0.34 

± 

0.09 

0.34 

± 

0.10 

3.9 

0.39 
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Table 3 — Continued 


ID 

Map ID 

R.A.(J2000) 

Decl.(J2000) 

gobs 

•^obs 

(mJy) 

ocorr 

0 1.2 mm 

SNR 

fsp 





(mJy) 



(1) 

(2) 




(3) 


(4) 

(5) 

(6) 

110 

60 

34.250927 

-4.802006 

0.37 

± 

0.11 

0.39 ± 0.13 

3.5 

0.68 

111 

61 

34.744495 

-4.855921 

0.35 

± 

0.10 

0.27 ± 0.12 

3.5 

0.70 

112 

61 

34.747295 

-4.858103 

0.37 

± 

0.10 

0.26 ± 0.11 

3.8 

0.45 

113 

62 

150.757401 

1.700786 

0.50 

± 

0.14 

0.57 ± 0.21 

3.5 

0.70 

114 

62 

150.758331 

1.703298 

0.42 

± 

0.12 

0.52 ± 0.17 

3.6 

0.60 

115 

62 

150.760971 

1.703856 

1.01 

± 

0.10 

1.17 ± 0.18 

9.7 

0.00 

116 

62 

150.762024 

1.698676 

0.57 

± 

0.16 

0.85 ± 0.24 

3.5 

0.69 

117 

62 

150.762161 

1.701653 

0.34 

± 

0.10 

0.37 ± 0.14 

3.5 

0.68 

118 

62 

150.762741 

1.700690 

0.43 

± 

0.12 

0.61 ± 0.17 

3.6 

0.62 

119 

62 

150.763138 

1.702965 

0.47 

± 

0.13 

0.70 ± 0.19 

3.6 

0.62 

120 

64 

34.587593 

-5.318766 

0.32 

± 

0.09 

0.21 ± 0.10 

3.6 

0.59 

121 

64 

34.589180 

-5.315989 

0.61 

± 

0.17 

0.69 ± 0.21 

3.5 

0.65 

122 

65 

34.387825 

-5.220457 

0.41 

± 

0.11 

0.42 ± 0.13 

3.7 

0.51 

123 

67 

197.860764 

-1.337225 

0.25 

± 

0.06 

0.03* ± 0.01 

4.3 

0.53 

124 

67 

197.862061 

-1.342226 

0.24 

± 

0.06 

0.00* ± 0.00 

4.3 

0.54 

125 

67 

197.864731 

-1.343368 

0.27 

± 

0.06 

0.03* ± 0.01 

4.7 

0.32 

126 

67 

197.869086 

-1.337143 

0.31 

± 

0.08 

0.08* ± 0.03 

4.1 

0.60 

127 

67 

197.871590 

-1.345783 

0.30 

± 

0.06 

0.02* ± 0.01 

5.1 

0.20 

128 

67 

197.874725 

-1.326138 

0.22 

± 

0.05 

0.02* ± 0.01 

4.3 

0.51 

129 

67 

197.881638 

-1.323127 

0.26 

± 

0.05 

0.05* ± 0.02 

4.8 

0.29 

130 

67 

197.882812 

-1.356609 

0.28 

± 

0.07 

0.05* ± 0.02 

4.1 

0.64 

131 

67 

197.883881 

-1.355241 

0.30 

± 

0.08 

0.06* ± 0.03 

4.0 

0.65 

132 

67 

197.884125 

-1.334885 

0.30 

± 

0.07 

0.13* ± 0.05 

4.2 

0.58 

133 

67 

197.886185 

-1.328756 

0.24 

± 

0.05 

0.04* ± 0.01 

5.1 

0.19 


S475 -7*625 > ~ 28 ^ r625 ~ ^ (^) 

where the source extraction and photometry are carried 
out with SExtractor. The external perturbation is cal¬ 
culated with the theoretical model under the assump tion 
that the perturbation is weak (e.g.. lKochancklli99li) . 

Using the posi tions of the multiple images prese nted 
in the literature (|Coe et al.ll 2010 l fPiego et' al.ll2015lh we 
optimize free parameters of the mass profiles based on 
the standard \ 2 minimization to determine the best-fit 
mass model. Figure [T| presents the best-fit mass model. 
The best-fit mass model achieves that all of the offsets 
between the model and observed positions of multiple 
images are within l'/0 in the image plane. We then cal¬ 
culate magnification factors p at the positions of our 
sources. For the p estimates, we assume our source red- 
shift of z = 2.5 that is the same as the one used in the 
flux scaling shown in Section [3] This assumption does 
not change our statistical results, because most of mm 
sources in the target fields at the high-galactic latitude 
should reside at z > 1 where the p values do not depend 
much on a redshift. Here we calculate the magnifica¬ 
tion factors for sources at z = 7.5 that correspond to 
the redshift of A1689-zDl, and evaluate the differences 
between the magnification factors for sources at z = 2.5 
(/H2.5) and z = 7.5 (/17.5). We find that the average dif¬ 
ference, Ap = (/i 7.5 — ^ 2 . 5 )/M 2 . 5 , is 0.36 ± 0.25 for our 11 
sources in A1689. Although the difference is not large, 
we include the lensing magnification difference of 0.36 
into the errors of the intrinsic flux estimates for our 11 
sources in A1689. Accordingly, we propagate these un¬ 
certainties to our major results such as number counts 
shown in Section 0] 

3.6. Survey Area 

We estimate survey areas of our data of A, B, and 
C, and present these estimates in Figure [5] The survey 


Note. — (1): Source ID. (2): Map ID that corresponds to 

the one in Table 2. (3) : Peak flux density of the SExtractor mea¬ 
surement (Section I3.lt with the primary beam correction at the 
observed wavelength. (4): The best-estimate source flux density 
at 1.2 mm. The source flux is estimated with the 2D Gaussian- 
fitting routine of imfit in CASA (Section 13.411 with the primary 
beam correction and the flux scaling to 1.2 mm (Section [3]). For 
the sources found in the cluster da ta, t he lensing magnification cor¬ 
rections are also applied (Section 13.511 . The lensing magnification 
factors are estimated at the ALMA flux peak positions. The error 
bar includes the random noise, 10% of ALMA system’s flux mea¬ 
surement ^uncertainty, and the lensing magnification uncertainty 
(Section mi ). (5): Signal-to-noise ratio of the peak flux density. 
(6): Spurious source rate that is defined by eq. [3] t Source iden¬ 
tified in the Band 7 map. Unless otherwise specified, the source is 
found in the Band 6 map. t The lensing magnification is corrected. 
The ID 124 source has 0.0020 ± 0.0010 mJy for the best-estimate 
source flux density at 1.2 mm with the primary-beam and lensing- 
magnification corrections. 

3.5. Mass Model 

We construct a mass model for the cluster data (C) of 
A1689 at z = 0.183. We make the mass model with 
the pa r ametr ic gravitational lensing pac kage GLAFIC 
( OerurU [20T0lf m the same manner as llshigaki ct al.1 
(12015 1. Our mass model consists of three types of 
mass distributions: cluster-scale halos, cluster member 
galaxy halos, and external perturbation. For evalu¬ 
ating the mass distributions, we make a galaxy cata¬ 
log of A1689, conducting source extractions in the op¬ 
tical images of 3475, ^25, and 1775 bands taken with 
Hubble Space Telescope (HST). The cluster-scale halos 
are estimated with three brightest member galaxies in 
the core of the cluster. The cluster member galaxies are 
selected by the color criteria of 

7*625 < 24, 

1 , 

3475 -7*625 < — To ( 7 ’625 _ 24) + 1.3, 
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Figure 5. Survey areas as a function of intrinsic flux densities. 
The dashed, dotted, and solid lines represent the survey areas of 
the field (A and B), cluster (C), and total (A, B, and C) datasets, 
respectively. 


areas are defined by the high sensitivity regions that are 
detailed in Section 13.11 Because the sensitivities of our 
ALMA maps are not spatially uniform, the survey areas 
depend on the flux densities. 

In the field data of A and B, the sensitivity of the pri¬ 
mary beam decreases with increasing radius from a map 
center. In other words, the detectable intrinsic flux den¬ 
sities (corrected for the primary beam attenuation) in¬ 
crease from the center to the edge of the data. First, we 
find the radius where sources with the intrinsic flux den¬ 
sities can be detected at our selection limit SNR. Then, 
we calculate areas with the radius given in each map, 
and sum up these areas to obtain the total survey areas. 

In the cluster data of C, the spatial distribution of the 
sensitivity is not expressed by a simple function of ra¬ 
dius, because the cluster data are taken by mosaic map¬ 
ping. Moreover, there are cluster-lensing magnification 
effects that allow us to detect intrinsically faint sources. 
We make an effective magnification map, multiplying, at 
each position, the sensitivity and the magnification factor 
estimated with the mass model of Section 13.51 We then 
calculate the survey area of the C data where a source 
with an intrinsic flux density can be observed above the 
selection limit. 

4. NUMBER COUNTS AND EBL 

4.1. Number Counts at 1.2 mm 

We derive differential number counts at 1.2 mm. For 
the number counts, we use the 122 sources identified in 
the ALMA Band 6 maps covering the 1.2 mm band, be¬ 
cause there is a possibility that the 11 sources in the Band 
7 would include unknown systematics in the flux scal¬ 
ing to 1.2 mm fluxes (Section 3). To derive the number 
counts, wc basically fo llow the methods used in the pre¬ 
vious ALMA studies fHatsukade et all 12013110 no et ahl 
l2014t iCarniani et al.ll2015fl . A contribution from an iden¬ 
tified source with an intrinsic flux density of S to the 


number counts, £, is determined by 


as) 


i-/ BP (g) 

c{s)A cB (sy 


(5) 


where C is the completeness, and A e g is the survey area. 
For C, A e ff, and / sp (eqU), we use the values derived in 
Section [3] Then, we calculate a sum of the contributions 
for each flux bin, 


n{S) 


ms) 

A log S” 


( 6 ) 


where A log S' is the scaling factor of 0.25 for the 1- 
dex width logarithmic differential number counts. We 
estimate the errors that include both Poisson statisti¬ 
cal errors of the source numbers and flux uncertainties 
of the identified sources. Because the numbers of our 
sources are s mall, we use th e Poisson uncertainty values 
presented in IGehrelsI (1198611 that are applicable for the 
small number statistics. The flux uncertainties are com¬ 
posed of the random noise and the ALMA system’s flux 
measurement uncertainty whose typical value is ~10% 
fSection [273]) . Moreover, for our 11 sources in the A1689 
data, the intrinsic flux estimates include the lensing mag¬ 
nification uncertainties (see Section 13.51) . To evaluate 
these flux uncertainties, we perform Monte-Carlo simu¬ 
lations. We make a mock catalog of the faint ALMA 
sources whose flux densities follow the Gaussian proba¬ 
bility distributions whose standard deviations are given 
by the combination of the random noise and the system 
(+ the lensing magnification) uncertainties. We obtain 
the number counts for each mock catalog in the same 
manner as those for our real sources. We repeat this 
processes 1000 times, and calculate the standard devi¬ 
ation of the number counts per flux bin caused by the 
flux uncertainties. Combining the standard deviation of 
the flux uncertainties with the Poisson errors, we finally 
obtain the lex uncertainties of the number counts. 

The differential number counts and the associated ler 
uncertainties are shown in Figure |6] and Table |4] With 
the technique same as those deriving the differential num¬ 
ber counts, we estimate the cumulative number counts 
that are summarized in Tabled 

In Figure [6] we find that the faintest data point at 
~ 0.002 mJy (red open circle) is composed of a source 
with a very high magnification factor of p > 100. The 
highly magnified sources would potentially include large 
systematic uncertainties originated from the mass mod¬ 
eling. Moreover, the faintest data point has the large 
error bar that does not provide important constraints on 
the number count measurements. Thus, we do not use 
the faintest data point in the following discussion. Ex¬ 
cept for this very high magnification factor object, the 
faintest and brightest sources in our sample have intrin¬ 
sic fluxes of 0.018 and 1.2 mJy, respectively. Thus, our 
study covers the flux density range of ~ 0.02 — 1 mJy. 


4.2. Comparison with Previous Number Count 
Measurements 

We compare our number counts to previous measure¬ 
ments at the ~ 1.2 mm band. In the flux range close to 
our study, there are two types of previous studies, blank- 
field observations with ALMA and cluster observations 
for lensed sources with single-dish telescopes. For the 
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Figure 6. Differential number counts at 1.2 mm. The red filled circles are our number counts derived from our faint ALMA sources. 
The red curve and the gray region denote the best-fit Schechter function and the associated lcr error. In the Schecht er function fitting, 
we use our number counts and the previous measur ements shown with the black filled stars (Scl2; I Scott et al .112013) and the magenta 
filled diamonds (Cal5 1.1 mm; ICarniani et all 120151 ). Note that the 1.1 mm results of Cal5 are obtained from ALMA sources, none of 
which are used in our number count measurements. Because the rest of the existing ALMA studies include some ALMA data covered by 
our stu dy, for our Schecht er function fitting we do not include the previous A LMA measurement s shown with the magenta open squares 
(Kal3; [Karim et al. 201 3]), op en inverse triangles (H13; Hatsukadc et al. 2013), open stars (014; R3no e t al. 20l |j), open diamonds (Cal5 
1.3 mm; [Carniani et al. 2015]), open triangles (Sil5; [Simpson et al. 2015a|), and open pentagons (Otl5: lOteo et all 120151) . Note that the 
faintest data point of our study (red open circle) is also removed from the sample of our Schether function fitting, due to the possible large 
systematic uncertainties. The pale-green region presents a flux range where no reliable number-count data bin exists (see text). We do not 
use the faintest data point of Scl2 (black open star) for our fitting, because this data point has unrealisticly small errors. The dashed and 
dotted curves show the model predictions for number counts based on co smological hydrodynamic simulations with GADGET-3 (Shl2; 
IShimizu et al.lf20121) and semi-analytic models (M15; IMakiva et al.l [20151) , respectively. These model predictions are roughly consistent 
with our results. For the previous measurements of the number counts at 1.1 mm and 1.3 mm different from our 1.2 mm band, we scale 
the flux densities with the methods shown in Section [3] Since the flux scaling factors for < 1 mm wavelengths would contain relatively 
large systematic uncertainties, we do not use the < 1 mm band data, for our Sche chter function fitting, that are presented with the open 
triangles (KnQ8: fKnudsen et al.l[2008fl and open hexagons (Chlfli fChen et al.H2013l ). 


measurements of 1.1 and 1.3 mm bands that are differ¬ 
ent from 1.2 mm band, we scale the flux densities with 
the methods described in Section [3] 

4.2.1. Blank-Field Observations 

To deriv e numb er counts. IHatsukade et al.1 (120131) and 
10no et ahl (120141) use 20 ALMA maps in one blank 
field of SXDS and 10 ALMA maps in 10 independent 
fields, respectively. ICarniani et al.1 ( 2015 ) obtain number 
counts with the 18 ALMA maps composed of 9 indepen¬ 
dent fields at 1.1 mm band and the 9 ALMA maps in 
one blank field of Cosmic Evolution Survey (COSMOS; 
iScoville et al.ll2007l ) at 1 .2 mm band. After the flux den¬ 


sity scaling, all of these previous number counts agree 
with our measurements within the lcr uncertainties. Be¬ 
cause the amount of our ALMA data is much larger than 
the previous studies, the statistical and systematic un¬ 
certainties of the number counts measurements are sig¬ 
nificantly smaller than those of previous studies. 

4.2.2. Cluster Observations for Lensed Sources 

Single-dish observations of SCUBA and SCUBA-2 for 
lensed sources behind galaxy clusters successfully reach 
the intrinsic flux limits compa rable to the p r evious 
ALMA blank field observations. iKnudsen ct al.1 (|2008D 
conduct SCUBA observations for 12 clusters, and de- 
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Table 4 

Differential and Cumulative Number 
Counts at 1.2 mm 


s 

l°g(n) 

(Vdiff 

(mJy) 

([A logS = lRMeg- 2 ) 


( 0 . 002 ) 

7 1 +0.5 

1 

0.020 

c c+0.3 
°-°-0.5 

2 

0.036 

C e+0.3 

O.D — 0.4 

6 

0.063 

c q+ 0.2 

O.0_o.4 

15 

0.112 

c 1 + 0.2 
°- 1 - 0.2 

25 

0.200 

4 9+5’i 
^• y - 0.2 

27 

0.356 

4 7+0-1 

1 - 0.2 

29 

0.633 

4 9+0.2 
^• z - 0.2 

15 

1.124 

3.4±°;| 

2 

S 

l°g(n(> S)) 

-/Vcum 

(mJy) 

(deg" 2 ) 


( 0 . 002 ) 

6 . 6^11 

122 

0.015 

k c+0.3 

O. 0 _o 4 

121 

0.027 

c q+ 0.2 

°-°-0.3 

119 

0.047 

c q+0.2 

O,U - 0.2 

113 

0.084 

^•°- 0.2 

98 

0.150 

4 6 +01 

73 

0.267 

4 2 “*" 0-2 
^• Z - 0.2 

46 

0.474 

q 7 + 0.2 

' —0.3 

17 

0.843 

9 0 +O .6 

Z -S-0.5 

2 


Note. — The lcr uncertainties are esti¬ 
mated from the combination of the number- 
count Poisson statistical errors and the flux 
uncertainties (see text). The column of JVdiff 
presents the numbers of our 1.2 mm sources in 
each flux bin, while the one of JV cum denotes 
the cumulative source numbers down to the 
flux densities. For our analysis, we do not use 
the faintest bin data indicated with the paren¬ 
theses (see text). 


rive the number counts at 850 fim with the sample of 
15 gravitationally lensed faint sources with flux densi¬ 
ties below the SCUBA’s blank - field c onfusion limit (~ 2 
mJy at 850 p m). I Chen et all (120131) observe 2 clusters 
as well as 3 blank fields with SC UBA-2, and obtain th e 
number counts similar to those of lKnudsen et al.1 ( 2008 1. 
The number counts from these previous studies are also 
shown in Figured These previous results generally agree 
with our results within the lcr errors. Although the num¬ 
ber counts of the previous studies are slightly higher than 
ours systematically at < 0.1 mJy, these small systematic 
differences would be explained by the uncertainties of the 
flux scaling from 850pm to 1.2 mm. 

4.3. Contributions to the EBL 

In Section 14.11 we obtain the 1.2-mm number counts 
with our faint ALMA source. Because our sample lacks 
sources brighter than 1.2 mJy, we use previous studies of 
bright sources in the literature to investigate the source 
contributions to the IR EBL. There are a number of stud¬ 
ies that investigate the bright sources (Figure [6]), but we 
use the number counts at 2 —10 mJy obtained from a sur¬ 
vey for six blank fields with a single -dish telescope instru¬ 
ment of AzTEC ( Scott et al.lf2012 ). which shows reliable 
number counts at 1.1 mm that is close to 1.2 mm band 
. Although there is another systematic survey for bright 
sources from the s ingle-dish observations with LABOCA 
dWeifi et al.H2009lb we do not use their data for our anal¬ 


ysis, due to their observation wavelength of 870 pm far 
from our 1.2 mm band. It should be noted that the 
num ber counts of bright sources are still under debate 
(e.g..lKarim et al1l2013t IChen et aI1l2013HSimpson et al.l 

I2015al) . because the recent high-resolution ALMA obser¬ 

vations reveal mu l tiple sources in a beam of single-dish 
observations dHodee et al.l 120131) . Nevertheless, our re¬ 
sults below do not depend on the bright-source number 
counts, because the contribution from the bright sources 
is not large (Section [[]). 

The source contributions to the IR EBL can be calcu¬ 
lated by integrating the number counts down to a flux 
density limit. To characterize the shape of t he number 
count s, we use the Schechter function form (|Schechterl 
I1976T ) in the same manner as the previous studies, 

mds = 0 . (£) «p (-!)<j(|.), <7) 

where </*, S'*, and a. are the normalization, characteristic 
flux density, and faint-end slope power-law index, respec¬ 
tively. The logarithmic form of n(S)d(logS) = (fr(S)dS 
is given by 

n{S) = (In 10)</* (Jpj exp jpj . (8) 


Using eq. ©, we conduct y 2 -fitting to the differen¬ 
tial number counts derived from our and the previous 
observations shown with the filled symbols in Figure [6] 
that are our measurements and the AzTEC (IScott et all 
|20T1) and ALMA 1.1 mm (|Carniani et al.ll2015[) results. 
Here we do not use the m e asurements of Hatsukade et al.1 
(12011) . IQno et, all (Ewiib IQteo et al l (1201 5lb and the 1.3 
mm results of iCarniani et al.1 (|2015l) that are presented 
with the open symbols of Figure El because these stud¬ 
ies use ALMA data covered by our study or the sam¬ 
ple that does not meet our selection criteria (Section 
& We confirm that the entire ALMA 1.1 mm data 
of lCarniani et all (|2015h are different from ours. Includ¬ 
ing the bright number counts at 2 — 10 mJy (see above), 
we perform Schechter function fitting with our number 
counts a nd those from th e 1.1-mm ob servations with 
AzTEC (IScott et al.l 1 201 ‘3l and ALMA (ICarniani et al.1 
I2015T) . We vary three Schechter parameters, and search 
for the best-fit parameter set of (</>*,S», a) that minimizes 
X 2 . The best-fit function and parameter set are presented 
in Figure [6] and Table [H respectively. Figure |6] and the 
X 2 /dof of 16.9/19 (Table E3) indicate that the number 
counts are well represented by the Schechter function 
down to ~ 0.02 mJy. 

We also carry out fitting of the double-power law 
(DPL) function in the same manner as the Schechter 
function fitting. The DPL function has four fr ee parame¬ 
ters o f </? PL , S* DPL , a DPL , and /3 DPL (see, e.g JOno et al.l 
120141) . The best-fit parameters are presented in Table [5j 


The x 2 /dof value for the DPL fitting is nearly unity 
(12.8/18), and the DPL also well represents the num¬ 
ber counts. Because the Schechter and DPL functions 
show similarly good fitting results, we use the Schechter 
function fitting results for simplicity in the following dis¬ 
cussion. 

With the best-fit Schechter parameter set, we calculate 
the integrated flux densities, f s S </(S) dS , down to the 
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Table 5 

Best-fit Parameters 


Schechter 

s. 


a 

X^/dof 

(mJy) 

(10 3 deg" 2 ) 



(1) 

(2) 

(3) 

(4) 

9 

0.16 

1 54 + 0 - 29 
i,o ^- 0.26 

9 1 9 + 0.07 
- 0.06 

16.6/19 

- DFE - 


<t>*^ 

a ufu puru 

X7 dof 

(mJy) 

(10 3 deg 2 ) 



(5) 

(6) 

(7) (8) 


4 43+u.cs 4 

0 26”*” U ' U5 

U.ZD_ 0> 07 

9 OQ+ 0 .U 5 9/-J+U.48 

z.oy_oo6 °- zu -i.n 

12.8/18 


Note. — (1)-(3) : Best-fit parameter set for the Schechter 
function. (4): \ 2 over the degree of freedom. (5)-(8): Best-fit 
parameter set for the DPL function. 


1 


SFR m [M 0 yr *] (z=2.5) 
10 100 



Figure 7. Fraction of resolved EBL at 1.2 mm as a function of 
the detection limit, Siimit* The red line and the gray shade indi¬ 
cate the best estimate from our study and the associated errors, 
respectively. The magenta, blue, and bla ck lines present the es¬ 
timates from th e previous ALMA stu dies ([Hatsukade et al .1 120131 : 
IQno et aH f2014: Carniani et al. 2 0151) . Here we re-calculate the 
resolved f ractions of t h e pre vious studies with the bright source 
counts of I Scott et al.l (120120 . T he blue shade den otes the IR 
EBL measurements from COBE dFixsen et al.l 119980 . The right 
axis represents the absolute values of the integrated flux densities, 
fg ” S(j)(S)dS, that correspond to the resolved EBL fraction. The 
top axis ticks IR star-formatio n rates SFRir, corresponding to the 
1.2 mm flux (|Kennicuttlll998t) for the case that a source resides at 
z = 2.5. 


flux limit of Siimit. Because the faintest bin of our num¬ 
ber count data covers down to 0.015 mJy (Figure 0, we 
choose Siimit = 0.015 mJy. We calculate the integrated 
flux density to be 22.9^1 Jy deg -2 at 1.2 mm. Fig¬ 
ure [7] shows the integrated flux density. We evaluate the 
fraction contributed to the t otal I R EB L meas urement 
at 1.2 mm (22/g 4 Jy deg -2 ; iFixsen et alJ 11998 ) that is 
obtained by the observations of the COBE Far Infrared 
Absolute Spectrophotometer (FIRAS), and find that the 
resolved fraction in our study is 104125%- This result 
suggests that our study fully resolves the 1.2 mm EBL 
into the individual sources within the uncertainties. 

Our result also indicates that the EBL contribution of 
sources with < 0.02 mJy are negligibly small. It would 


suggest that there is 1) a flattening of number counts 
slope a, or 2) a truncation of the number counts at a flux 
«5trunc below ~ 0.02 mJy. Using our resolved fraction of 
/res = 104 / 25 %) we constrain a and Strunc i n the scenar¬ 
ios of 1) and 2). For the scenario of 1), we extrapolate 
the Schechter function from S'limit to 0 mJy with a power 
law of n(S ) oc sJ a °- s limit + 1), where ao_s Umit is a power 
law slope at 0 mJy to S’limit, and integrate the Schechter 
function with the power law extrapolation. We constrain 
the ao-Snmit value for the condition that the integrated 
flux density does not exceed the ler upper limit of the 
EBL measurement. We thus determine that the lower 
limit of ao-Snmit ^ —1.0. For the scenario of 2), we ex¬ 
trapolate our best-fit Schechter function below Siimit- We 
find Strun c — 0-004 mJy (c orresponding to SFJ?ir ~ 1.6 
M 0 yr -1 ; iKennicuttl 1 199811 where the ler lower error of 
/res reaches 100% of the EBL contribution fraction. Be¬ 
cause these 1) and 2) scenarios are two extreme cases, our 
results would suggest that the faint-end slope at < 0.02 
mJy is flatter than the best-fit Schechter function slope 
of a = —2.12/g;Qg (TableO, and that there is a number- 
counts truncation, if any, at < 0.02 mJy. 


5. CLUSTERING ANALYSIS 


5.1. Galaxy Bias of our Faint ALMA Sources 

We estimate a galaxy bias of our fai nt ALMA sources 
by the counts-in-cells technique (e.g. lAdelberger et alJ 
119981 ; Robertsonil2010H 


N 2 ay(z) ’ 


where c+ is the standard deviation of the detected source 
counts per field, ay is the matter variance averaged over 
the survey volumes V, and N is the average source counts 
per field. 

In this counts-in-cells analysis, we use 66 maps of the 
held data (A) and (B). Because the cluster data (C) in¬ 
cludes lensing effects of magnification and survey vol¬ 
ume distorti on a s well as the different spurious source 
rate (Section [372]) , we do not include the cluster data (C) 
in the sample of our counts-in-cells analysis. The sur¬ 
vey volume is calculate d by using the average effective 
area per held fSection 13.611 . For the redshift distribu¬ 
tion to calculate the volume, we again assume a median 
redshift of z — 2.5 (Section [3]) and the top hat func¬ 
tion covering z = 1 — 4. For a% and N, we do not use 
the simple detection numbers, but the corrected detec¬ 
tion numbers obtained with equation ([oj , which allows us 
to remove the effects of spurious sources. We calculate 
o' 2 - with the analy tic structure formation model (e.g., 
IMo fc White! 120(B ). To estimate the ler uncertainty of 
b g , we carry out Monte-Carlo simulations. We hrst make 
1000 mock catalogs of faint ALMA sources consisting of 
source counts that follow the Poission probability dis¬ 
tribution function whose average values agree with our 
observed source counts. Then, we calculate galaxy bias 
values with each mock catalog, and repeat this process 
for the 1000 mock catalogs. We define a ler error of our 
observational b g estimate by the 68 percentile of the b g 
distribution obtained from the 1000 mock catalog results. 
Because the ler error is large, we cannot obtain the mea¬ 
surement of b g but place a ler upper limit of b g < 3.5 
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Figure 8. Galaxy bias of our faint ALMA sources and the 
other high-2 galaxies. The red arrow shows the lcr up¬ 
per limit of our faint ALMA sources. The magenta arrows 
are the measurements and the upper limit for br ight SMGs, 
which a re calculate d from t he clusterin g resul ts of IWebb et al 
(2003); Blair^t_alJ (120041); IWeifi et al.l (120091) ; I Williams et al. 
(2011); IHickox et al.l (12012) . The orange and green squares de¬ 
note t he if-selected galaxies of DR Gs/pBzK and sBzK, respec¬ 
tively dGrazian et al.ll20061;IHavashi et al .1 1200 71; I Qua dri et'al.H20071 ; 
IBlanc et alJ 120081 : IFurusawa et alJ 120111 : I Lin et al.l 12012 ). The 

blue and cyan diamonds represent the UV-selected galaxies of 
LBGs (including BX/BM) and LAEs, respectiv e ly dOuchi et al.l 
200J^JAdelberg^et ajj[2QQ5| JOverzier et al .1120061 ; ILee et al.l 120061 ; 

Gawiser et alil2007l : IQuchi et al.l 120101 ). The solid curves indicate 

the galaxy biases of the dark halos with masses of 10 9 , 10 10 , 10 11 , 
10 12 , an d 10 1 3 Mr. . from bottom to top, which are predicted by the 
model of IShcth fe Tormcnl (119991 1 under the assumption of the one- 
to-one relation between galaxies and dark halos. The dotted curves 
present evolutionary tracks of dark halos for the galaxy conserving 
case, which assume that the m otion of ga laxy is purely caused by 
gravity with no galaxy merges l !Frvill996i ). 


for our faint ALMA sources. This lcr upper limit agrees 
with the previo us b n estimate for faint ALMA sources, 
b g < 4, given bv lOno et al.l (|2014f ). 

We estimate the dark halo masses of our faint ALMA 
sources with our result of b g < 3.5 with the a n alytic s 
structure formation model of iShcth fc Tormcnl ( 1999 ). 
assuming the one-to-one correspondence between galax¬ 
ies and dark halos. We obtain the lcr upper limit dark 
halo mass of < 5 x 10 12 M 0 . 

5.2. Comparison with Other Populations 

Figure [5] shows the galaxy biases of our faint ALMA 
sources and a variety of high-z galaxy populations; bright 
SMGs, Af-selected galaxies including passively evolving 
BzK galaxies (pBzKs) and star-forming BzK galaxies 
(sBzKs), distant red galaxies (DRGs), and UV-selected 
galaxies composed of BX/BM galaxies, LBGs, and LAEs. 
Some of these high-z galaxy studies shows no galaxy bias 
estimates, but only the galaxy correlation lengths ro, and 
the power-law indices of the two point spatial correlation 
function 7 . For these results, we calculate the galaxy 


biases in the same manner as lOno et akl (|2014f) from 
b g = erg g a i/<T8(z), where cr 8 (z) is a matter fluctuation 
in spheres of comoving radius of 8 h _1 Mpc, and (78 lga i is 
a galaxy fluctu ation. The va lue of cr gjga i is derived with 
(see eq. 7.72 in lPccblcslll993l l 

2 _ 72 ( r ° 

^8, gal - (3 _ 7 )( 4 _ 7 )( 6 _ 7 )27 ^8 h - 1 Mpc 

In Figure HI the galaxy bias estimates for the most 
of SMGs, DRGs, and pBzKs are higher than the up¬ 
per limit of our faint ALMA sources. On the other 
hand, AT-selected star-forming galaxies of sBzKs have the 
halo masses comparable with those of our faint ALMA 
sources. The UV-selected LBGs and LAEs also have the 
galaxy bias values similar to our faint ALMA sources. 
Thus, we find that our faint ALMA sources are not simi¬ 
lar to SMGs, DRGs, and pBzK with masses of > 1O 13 M 0 
but sBzKs, LBGs, and LAEs with masses < 1O 13 M 0 
(Figure 0). These clustering results suggest a possi¬ 
bility in a statistical sense that a large fraction of our 
faint ALMA sources could be mm counterparts of sBzKs, 
LBGs, and LAEs. 

6. MULTI-WAVELENGTH PROPERTIES OF THE FAINT 
ALMA SOURCES 

6 .1. Optical-NIR Counterparts 

We investigate optical-NIR counterparts of our faint 
ALMA sources, and characterize the sources with the 
multi-wavelength properties, especially the optical-NIR 
colors. We search for the optical-NIR counterparts in 
the regions of SXDS and A1689 that have rich multi¬ 
wavelength data. In these two regions, there are a total of 
65 faint ALMA sources; 54 and 11 sources in SXDS and 
A1689 regions, respectively. For the optical-NIR counter¬ 
part studies so far conducted, this is the largest sample of 
faint mm sources down to ~ 0.02 mJy, while there exist 
studies, such as ALESS and COSMOS, tha t use a com¬ 
parably large but relative ly bright sample (|Hodge et akl 
Effl IScoville et al.lEfTl . 

Table [6] summarizes the multi-wavelength data of 
SXDS and A1689 used in this study. The multi¬ 
wavelength data cover from X-ray to far-infrared (FIR) 
radiation. In the SXDS region, we use published pho- 
tometry catalogs in 0.5-10 keV (lUeda et al.l I2008T) . u* 
Foucaud et al. in preparation), B, V, R, i', z' 
Fnrnsawa et hklFhliul: . J, H, K (UKIDSS DR10; Al- 
maini et al. in preparation), 3.6 — 24 /an (Spitzer 
UKIDSS Ultra Deep Survey; SpUDS), 250 — 500 /im 
(Herschel Multi-ti ered Extragalac tic Survey; HerMES , 
lOliver et akl 120121) . and 1.4 GHz dSimnson et al.l [2006h 
as well as the phot ometric redshift (z p hot) catalog 
(jWilliams et al.l 12009 1. In the A1689 region, we in¬ 
vestigate published photometry catalogs of <7475, ^25, 
* 775 ; ^ 850 ; <7no, Hi 60 (Hubble Source Catalog; HSC), 
3.6 — 24 /rm ( Spitzer Enhanced Imaging Products; 
SEIP), 250 — 500 /im (HerMES), and the photomet¬ 
ric re dshift z n hntn and spectroscopic redsh i ft z sr>Rr cat¬ 
alogs (jLimousin et al.ll2007i : iCoe et al.ll20Tcil : iDiego et akl 
1201511 . 

We use the u* to K ( g 1 to H) band images in the SXDS 
(A1689) region for the optical-NIR (0.4 — 2 /im) coun¬ 
terpart identification, because these data have a good 
spatial resolution. The optical-NIR counterparts are de- 
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fined by the criterion that the distance between optical- 
NIR and ALMA source centers is within a l'/O. If an 
ALMA source meets this criterion in any of optical-NIR 
bands, we regard that the ALMA source has an optical- 
NIR counterpart. The l'/O radius is chosen, because the 
absolute positional accuracy (APA) of ~ l'/O is guaran¬ 
teed in our A LMA data (cf. op tical-NIR. data for < 0'.'3; 
see, e.g.. iFurusawa et aLll2008lh Regarding the ALMA’s 
APA, the knowledge base page of ALMA Science portal 
presents that the APA is smaller than the synthesized 
beam widtl0. Because the synthesized beam widths of 
our data are < l'/O, the l'/O radius is applied. 

Using the SXDS+A1689 source catalog of the 65 faint 
ALMA sources, we identify a total of 17 optical-NIR 
counterparts. Fifteen and two sources are identified in 
SXDS and A1689 and referred to as Sl-15 and Al-2, re¬ 
spectively. Figure fl5l shows the postage stamps of these 
identified sources, and Table Q] summarizes their photo¬ 
metric properties. The photometry values in Table [7] are 
total magnitudes. The u* magnitudes are MAG_AUT0 val¬ 
ues given by our SExtractor photometry. Similarly, 
BVRri'z' magn i tudes are MAG_AUT0 values listed in the 
IFurusawa et al. (2008)’s catalogs. The Spitzer photome¬ 
try values are estimated from the SEIP aperture flux den¬ 
sities with the aperture correction Q while the HST pho¬ 
tometry values (for Al and A2) are the total magnitudes 
calculated from MAG_APER values of our SExtractor 
measurements and the aperture correction @. Note that 
there is an offset of l'/2 (l'/06) between the flux peaks of 
the optical-NIR counterpart and S9 (S15) ALMA source. 
Because S9 and S15 have an extended optical-NIR profile 
whose major component falls in the ALMA emitting re¬ 
gion, the optical-NIR source is selected as a counterpart 
of the ALMA source. 

Because spurious sources are included in the 65 faint 
ALMA sources, we calculate the real ALMA source num¬ 
ber in the SXDS (A1689) region. In the SXDS (A1689) 
region, we subtract the expected spurious-source num¬ 
ber of 31 (5) from the SXDS (A1689) source number 
of 54 (11), and estimate the real ALMA source num¬ 
ber to be 23 (6). The total real ALMA source num¬ 
ber in these two regions is 29 (= 23 + 6). Thus, the 
counterpart-identification fraction that is the ratio of the 
total optical-NIR counterpart ALMA sources to the real 
ALMA sources is estimated to be 59% (=17/29). We 
estimate the probability of the chance projection that 
a foreground or background object is located by chance 
in the l'/O-radius circl e of an ALMA sourc e, and obtain 
~ 10% in the P-value dDownes et al.lfl986lh Multiplying 
this probability with the number of optical-NIR coun¬ 
terparts, 17, we expect that ~ 2 objects would be the 
foreground/background contamination. Even if we in¬ 
clude this contamination effect, a half of our faint ALMA 
sources, 52% (= [17 — 2]/29), have an optical-NIR coun¬ 
terpart. 

Interestingly, there is a difference of the counterpart- 
identification fractions between SXDS and A1689 re- 

6 The APA is larger than 0'.'1 that is originated from the un¬ 
certainties of phase calibrator positions and the baseline lengths of 
antennas. 

http://irsa.ipac.caltech.edu/data/SPITZER/docs/irac/ 
iracinstrumenthandbook / 29 / 

http://archive.stsci.edu/hst/hsc/help/FAQ/aperture 
.corrections. txt 


gions, although the depths of the optical-NIR data are 
comparable in these two regions (27 — 28 mag in op¬ 
tical and 25 — 26 mag in NIR; see Table [6J. The 
counterpart-identification fraction in the SXDS region 
is 65% (= 15/23), while the one in the A1689 region 
is only 33% (= 2/6). Including the ler statistical un¬ 
certainties, these fractions are 65 ± 22% and 33 ± 30% 
for the SXDS and A1689 regions, respectively. We care¬ 
fully examine whether this difference is just an artifact 
originated from, for examples, a systematically high spu¬ 
rious rate and cluster member galaxy obscurations in the 
cluster region. However, t here are no hi nts o f notable ar¬ 
tificial effects in the data. iChen et al.l (|2014f ) also report 
a similarly low-identification rate of MIR/radio counter¬ 
parts in A1689 based on their observations with SMA. 
It should be noted that the intrinsic mm flux densities 
of our sources in the A1689 region are ~ 0.01 — 0.1 mJy 
(magnification corrected) that is about an order of mag¬ 
nitude fainter than those of our sources in the SXDS 
region, ~ 0.1 — 1 mJy. 

There are two possible explanations for the difference 
between SXDS and A1689. The first explanation is made 
by the negative k-correction. Due to the faint flux densi¬ 
ties, the mm sources of the A1689 would reside at a red- 
shift higher than the SXDS sources on average, since the 
abundance of the dusty population decreases towards the 
epoch of the first-generation galaxies. Because the effect 
of the negative k-correction is significant, the ratio of mm 
to optical-NIR fluxes is higher for high-z sources than 
low -z sources. In this case, the counterpart-identification 
fraction would be small in the A1689 region. The sec¬ 
ond explanation is the A1689’s lensing distortion effects 
in conjunction with the intrinsic offsets betw een mm and 
optica l-NIR. emitting regions within a galaxy. IChen et all 
( 2015 1 report that there is an offset between submm and 
optical-NIR emitting regions typically by ~ O'.'4. If there 
is an intrinsic ~ 0'/4 offset in the source plane of A1689, 
the image distortion effects of lensing would magnify the 
offset to > l'/O. We thus estimate the probability that 
the offsets between mm and optical-NIR sources in the 
source plane is magnified to > l'/O in the image plane by 
Monte-Carlo simulations. We place 100 artificial sources 
in the source plane of A1689 at around an ALMA source 
in A1689. The positions of these artificial sources fol¬ 
low the Gaussian distributions whose standard deviation 
is O'.'4. Using our mass model, we calculate the offsets 
between the ALMA source and the artificial sources in 
the image plane. We repeat this process for 11 ALMA 
sources in A1689. We find that ~ 50% of the artificial 
sources show > l'/O offsets from the ALMA source in the 
image plane. Because the number of optical counterparts 
in A1689 is 2, our simulation results suggest that a to¬ 
tal of 4 (= 2/[50%]) optical counterparts should exist in 
the case of no lensing distortion effects (+the intrinsic 
offsets). Thus, the true counterpart-identification frac¬ 
tion of A1689 is 67%(= 4/6) that is very similar to the 
counterpart-identification fraction of SXDS (65 ± 22%). 
If the lensing distortion effects (+the intrinsic offsets) 
would provide 2 (=4 — 2) sources whose optical counter¬ 
parts have an offset from the ALMA sources by > l'/O 
in the image plane, we would find such ALMA sources 
whose optical counterparts are located at the position 
just beyond l'/O in our data. We refer the distribution 
of the artificial optical counterparts from our Simula- 
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tions, and search for potential optical-NIR counterparts 
whose distances from the ALMA sources are > l'.'O. We 
find such potential optical counterparts in 4 sources of 
A1689, some of which may be real optical counterparts, 
and confirm the hypothesis that the combination of the 
lensing distortion effects and the intrinsic offsets can be 
the major reason for the difference of the counterpart- 
identification fractions between the cluster of A1689 and 
the blank field of SXDS. 

Table[7]lists the photometry of the Spitzer (3.6—24 pm) 
and Herschel (250 — 500 pm) data. Because the spatial 
resolutions of the Spitzer and Herschel data are poor, 
we determine the Spitzer and Herschel counterparts by 
visual inspection that allows sources residing at a po¬ 
sition slightly beyond the 1" search radius. Although 
no sources are identified in the Herschel bands, about 
a half of our faint ALMA sources are detected in the 
Spitzer data. These results would indicate that our faint 
ALMA sources are very high redshift sources (z > 4). 
However, it should be noted that the Herschel data can 
detect dusty starbursts whose SFRs are very high, > a 
few lOO-Mgyr -1 even at z ~ 2 — 3. 

Here we examine whether we simply miss red spectral 
energy distribution (SED) objects that are not detected 
in optical-NIR (0.4 — 2 pm) bands but only in Spitzer 
and Herschel bands. We search for Spitzer and Herschel 
counterparts of our faint ALMA sources with no optical- 
NIR counterparts, adopting the visual inspection tech¬ 
nique described above. We find no Spitzer and Herschel 
counterparts that have undetectable optical-NIR fluxes. 
Again, these results do not suggest that our faint ALMA 
sources are dusty starbursts. 

We cross-match the catalogs of our faint ALMA 
sources and the radio 1.4 GHz and X-ray 0.5-10 keV 
sources. We find no counterparts within a l'/0 — 3'/0 
radius around our faint ALMA sources. Thus, our faint 
ALMA sources have no signatures of AGN. 

To summarize these counterpart identification results 
with the Spitzer and Herschel as well as radio and X- 
ray data, there are signs of neither dusty starbursts nor 
AGNs in our faint ALMA sources. Note that these cross¬ 
match analysis is completed only in SXDS because there 
are no radio and X-ray data available in A1689. 

6 .2. Lensed ALMA Sources in A1689 

We identify two optical counterparts of Al and A2 that 
are lensed by the A1689 cluster. The optical image at the 
position of Al is presented in Figure [2] There is a small 
positional offset of ~ 0'/8 between the optical and mm 
emission, but this small difference can be explained by 
the positional uncertainty of the ALMA image (Section 

EDB 

This optical counterpart of Al is one of three known 
multiple images that is dubbed 5.2, while t he rest of two 
multi ple images are called 5.1 and 5.3 (|Limousin et all 
[2007 1. We investigate our ALMA image at the optical 
positions of 5.1 and 5.3 sources. In the ALMA image, the 
5.1 source is marginally detected at the 3er level, while 
the 5.3 source has no signals above the 3cr level. Table 

9 There remains a possibility that the offset would be a real 
positional difference between the optical and mm emitting regions. 
However, it requires ALMA data analysis for a positional accuracy 
better than this study. 


[8] summarizes the physical properties of these three mul¬ 
tiple images. The difference of mm emissivities between 
these multiple images is discussed in Section RTT1 

Although the 5.2 source, i.e. Al, does not have a spec¬ 
troscopic redshift, the redshift of the 5.1 and 5.3 sources 
is z = 2.60 deter mined b y spectroscopy (Richard et al. in 
preparation, see iLimousin et ah! 120071) . We thus regard 
that the redshift of Al is z = 2.60. We calculate the 
intrinsic 1.2 mm flux of Al from our 1.3 mm data with 
the flux scaling factor from 1.3 mm to 1.2 mm (Section 
[3]), and obtain 0.020 ± 0.008 mJy that is corrected for 
the gravitational lensing mag nification of p = 22.8 and 
the primary beam sensitivity. iChen et all (120141) report 
that the 5.1 and 5.3 sources are not detected above the 
3cr level by their Submillimcter Arr ay (SMA ) 870 (im 
observations, but only the 5.2 source. [Chen et ah 1 (120141 ) 
measure the flux of the 5.2 source, and find the intrinsic 
flux of 5.2 source is 0.085 ±0.035 mJy at 870 /im, which 
corresponds to 0.033 ± 0.014 mJy at 1.2 mm with the 
fc-correction from 870 /jrn to 1.2 mm (Section [3]). This 
1.2 mm flux of 0.033 ± 0.014 mJy is consistent with our 
measurement of 0.020 ± 0.008 mJy within the ler uncer¬ 
tainties. 

The optical-NIR counterpart of A2 is identified in 
the optical bands with a lensing magnification factor of 
p = 4.9. The intrinsic 1.2 mm flux of A2 is estimated 
to be 0.077 ± 0.035 mJy that are corrected for the gravi¬ 
tational lensing magnification and the primary beam at¬ 
tenuation. Neither z spe c nor z p hoto of A2 are available in 
the literature. 

We compare these A1689 sources with optical- 
NIR counterparts (Table |SJ) and the p reviously - known 
lense d sources, Cosmic Evla sh fe.g.. iSwinbank et al.l 
l2010h and SMMJ16358 (e.g., iKneib et all l2004h . iden- 
tified by single-dish observations. Cosmic Eyelash 
and SMMJ16358 have intrinsic SFRs of ~ 210 and 
500M e yr _1 , respectively. These SFRs correspond to 0.5 
and 1.3 mJy at 1.2 mm under the as sumption of the 
SFR — Lir, relation dKcnnicuttl I1998D . On the other 
hand, our optical-NIR counterparts have intrinsic fluxes 
of ~ 0.02 — 0.08 mJy at 1.2 mm. Thus, our ALMA study 
has identified sources with intrinsic fluxes > 10 times 
fainter than these lensed sources. 

Although there are no optical-NIR counterparts, we 
find a very strongly lensed source candidate, A3, that 
shows ALMA flux double peaks dubbed A3a and A3b. 
The upper panel of Figure [TO] shows the ALMA image 
of A3. The A3 source is very close to the critical line 
for a ^ = 3 source /Figure flOl) . and the critical line is 
located between the double peaks of A3. More quanti¬ 
tatively, we investigate whether these double peaks can 
be explained by multiple images of a single lensed source 
behind A1689. Using our mass model, we calculate po¬ 
sitions of multiple images made by a lensed source at 
z = 2 — 6 (cyan lines in the upper panel of Figure [1(11). 
We find that multiple images at 2 ~ 3 reproduce the 
positions of the double peaks very well. The best-fit red¬ 
shift is z = 3.07, and the magnification factor is very 
high, 159. The bottom panel of Figure |Tt)] is the HST 
optical images. No optical counterparts of A3 are found. 
We summarize the properties of A3 in Table [9] 
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Table 6 

Multi-wavelength Data 


Region (Ncounter) Telescope/Instrument 

Band (5 a limit in mag) 

Reference 

(i) 


(2) 

(3) 

SXDS (15) 

XMM-Newton/EPICs, pn, MOS 0.5-10 keV (6 — 50 x 10 16 erg 

cm -2 s _1 ) (a) 


CFHT/MegaCam 

u*( 27.7) 

(b) 


Subaru/Suprime-Cam 

B (28.4) 

(c) 


Subaru/Suprime-Cam 

V (28.0) 

(c) 


Subaru/Suprime-Cam 

He (27.8) 

(c) 


Subaru/Suprime-Cam 

i' (27.7) 

(c) 


Subaru/Suprime-Cam 

z' (26.7) 

(c) 


UKIRT/WFCAM 

J (25.6) 

(d) 


UKIRT/WFCAM 

H (24.8) 

(d) 


UKIRT/WFCAM 

K (25.4) 

(d) 


Spitzer/ IRAC 

3.6pm (24.7) 

(e) 


Spitzer/ IRAC 

4.5pm (24.6) 

(e) 


Spitzer/IRAC 

5.8pm (22.3) 

(f) 


Spitzer/ IRAC 

8.0/rni (22.3) 

(0 


Spitzer/ MIPS 

24pm (18.0) 

(0 


Herschel/ SPIRE 

250pm (13.8) 

(g) 


Herschel/ SPIRE 

350pm (14.0) 

(g) 


Herschel/ SPIRE 

500pm (13.6) 

(g) 


VLA 

1.4 GHz (20.0) 

(h) 

A1689 (2) 

HST/ACS 

9475 (28.6) 

1 


HST/ACS 

7625 (28.3) 

(i) 


HST/ACS 

*775 (28.2) 

(i) 


HST/ACS 

2850 (27.8) 

(i) 


HST/mCMOS 

Aio (26.2) 

(j) 


HST/ WFC3 

H±6o (26.3 - 28.1) 

(k) 


Spitzer/ IRAC 

3.6pm (24.8) 

(1) 


Spitzer/YRAC 

4.5pm (24.3) 

(1) 


Spitzer/ IRAC 

5.8pm (23.4) 

(1) 


Spitzer/IRAC 

8.0/im (22.9) 

(1) 


Spitzer/ MIPS 

24pm (19.5) 

(1) 


Herschel/SPIRE 

250pm (14.0) 

(g) 


Herschel/ SPIRE 

350pm (14.2) 

(g) 


Herschel/ SPIRE 

500pm (13.8) 

(g) 

Note. - 

Columns: (1) Name of the region. 

The number, lV r onnter, in the parenseses 

indicates the number of 

the optical-NIR counterparts. (2) Limiting magnitude defined by 5cr sky noise in a circular aperture with a radius 

of PSF FWHM. (3) Reference or source catalogue: fa) lUeda et al.l (120080: (b) Foucaud et 

al. (in preparation); (c) 

Furusawa et al. (20081): (d) UKIDSS fO. Almaini 

et al. in preparation); (e) SEDS dAshbv 

et al. 201311: (l) SdUDS: 

(s) HerMES 

( Oliver et al.1120121): fh) Simpson et 

al.l (1200GT): (i) Alavi et al. (2014); fi) IB< 

auwens et al.l (20091): fk) 


iBouwens et all 


a http://irsa.ipac.caltech.edu/data/SPITZER/Enhanced/SEIP/overview.html 


6.3. Color and Luminosity Properties of the 
Optical-NIR Counterparts 

We investigate color (i.e. SED) and luminosity prop¬ 
erties of the optical-NIR counterparts, Sl-15 and Al- 
2, to understand the population of the faint ALMA 
sources with the detectable optical-NIR continuum emis¬ 
sion. First of all, we refer the c atalog o f p hoto-est i- 
mated from the optical-NIR SEDs (jWilliams et al.ll2009ft . 
and find that Sll, S12, S15 have photo -2 values of 
1 .54+QQ2, 1.57 ^q;o2j an d 1.45lo;g2, respectively. How¬ 
ever, the photo -2 values of the S4 and S6 sources are 
0 .771 o'q3 and 0.67 /qq 2, respectively. We thus regard 
that the S4 and S6 sources reside at low redshift (z < 1), 
and find that ~ 10% (= 2/17) optical-NIR counterparts 
are low -2 objects. Because the majority of mm sources 
are th ought to be high -2 galaxies and AGNs (|Casev et al.l 
12014ft . the rest of 15 (= 17 — 2) optical-NIR counterparts 
are candidates of high -2 galaxies. We thus compare col¬ 
ors and luminosities of these high -2 galaxy candidates 
with those of high -2 sources. The major high -2 pop¬ 
ulations are SMGs, LBGs (including BX/BM), DRGs, 
BzKs, LAEs, and AGNs. Since the SMGs, DRGs, and 
LAEs generally have weak optical continua that cannot 


be reliably compared with our optical-NIR counterparts, 
we focus on the comparisons with LBGs (BX/BM) and 
BzKs. 


6.3.1. Comparison with the LBG BX/BM Populations 

The BX and BM galaxies, generally included in the 
LBG population, reside in the redshift ra nges of 2.0 < 
2 < 2.5 and 1.5 < 2 < 2.0, respectively (jSteidel et al.l 
12004ft . Similarly, there are LBGs at 2 ~ 3. These galax¬ 
ies are defined by the selection criteria of U n — G and 
G — R colors. We apply these color selection criteria to 
our optical-NIR counterparts, but the photometry bands 
of our data are different fr om the (7 » , G a nd R bands. 
Following the procedures of iLv et ah ' (1 2011 7) , we alterna¬ 
tively use a color set of U — BV and BV — R c i'. The BV 
and R c i' photometry are defined by 


BV = — 2.51og 
R c i' = —2.51og 


XlfB + (1 - Xify) 
3630 pJy 
X 2 fa + (1 - X 2 fi>) 
3630 [iJy 


( 11 ) 


where fx is the flux density per unit frequency (erg s 1 
cm -2 Hz -1 ) in a band ’X’ corresponding to B , V, R c , 












































Table 7 

Photometry of the Optical-NIR Counterparts 


SXDS 


ID 

u* 

B 

V 

Rc 

i' 

z' 

J 

H 

K 

3.6pm 

4.5pm 

5.8pm 

8.0pm 

24pm 

250pm 

350/im 

500/im 

SI 


24.46 

24.24 

23.85 

23.52 

22.88 




20.02 

19.94 

20.11 

20.12 

> 18.56 

> 

14.35 

> 

14.55 

> 14.15 

S2 


25.91 

26.18 

26.59 

26.71 

25.01 




> 25.27 

> 25.16 

> 22.85 

> 22.81 

> 18.56 

> 

14.35 

> 

14.55 

> 14.15 

S3 

27.24 

26.69 

26.35 

25.63 

25.29 

25.25 

24.80 

24.43 

24.41 

> 25.27 

> 25.16 

> 22.85 

> 22.81 

> 18.56 

> 

14.35 

> 

14.55 

> 14.15 

S4 

25.76 

25.61 

25.30 

25.00 

24.47 

24.05 

23.21 

22.92 

22.66 

21.68 

21.84 

> 22.85 

> 22.81 

> 18.56 

> 

14.35 

> 

14.55 

> 14.15 

S5 

26.31 

25.11 

24.70 

24.76 

24.58 

24.75 

24.53 

24.05 

24.22 

> 25.27 

> 25.16 

> 22.85 

> 22.81 

> 18.56 

> 

14.35 

> 

14.55 

> 14.15 

S6 

25.69 

24.84 

23.76 

22.79 

21.96 

21.55 

20.81 

20.42 

20.08 

19.61 

20.06 

19.95 

20.24 

> 18.56 

> 

14.35 

> 

14.55 

> 14.15 

S7 

> 28.23 

27.14 

27.42 

27.30 

27.32 

26.62 

> 26.18 

> 25.31 

> 25.94 

> 25.27 

> 25.16 

> 22.85 

> 22.81 

> 18.56 

> 

14.35 

> 

14.55 

> 14.15 

S8t 

24.57 

> 28.93 

> 28.56 

> 28.33 

> 28.24 

> 27.22 

> 26.18 

> 25.31 

> 25.94 

> 25.27 

> 25.16 

> 22.85 

> 22.81 

> 18.56 

> 

14.35 

> 

14.55 

> 14.15 

S9tt 

> 28.23 

25.30 

23.86 

22.96 

21.63 

20.93 

20.24 

20.22 

20.45 

20.40 

20.78 

20.46 

> 22.81 

> 18.56 

> 

14.35 

> 

14.55 

> 14.15 

S10 

> 28.23 

26.68 

26.25 

26.25 

25.65 

25.07 

25.90 

24.60 

25.67 

> 25.27 

> 25.16 

> 22.85 

> 22.81 

> 18.56 

> 

14.35 

> 

14.55 

> 14.15 

Sll 

26.00 

25.14 

24.66 

24.53 

24.02 

23.62 

22.36 

21.76 

21.26 

19.91 

19.69 

19.31 

19.51 

15.95 

> 

14.35 

> 

14.55 

> 14.15 

S12 

23.45 

23.22 

22.96 

22.94 

22.65 

22.45 

21.55 

21.18 

20.88 

19.89 

19.66 

19.35 

19.62 

16.66 

> 

14.35 

> 

14.55 

> 14.15 

S13 

> 28.23 

> 28.93 

28.50 

> 28.33 

> 28.24 

> 27.22 

> 26.18 

> 25.31 

25.22 

22.55 

22.73 

22.33 

> 22.81 

> 18.56 

> 

14.35 

> 

14.55 

> 14.15 

S14 

> 28.23 

25.96 

26.09 

25.68 

> 28.24 

> 27.22 

> 26.18 

> 25.31 

> 25.94 

> 25.27 

> 25.16 

> 22.85 

> 22.81 

> 18.56 

> 

14.35 

> 

14.55 

> 14.15 

S15 

23.87 

23.49 

23.40 

23.08 

22.84 

22.42 

21.72 

21.25 

21.16 

20.27 

20.12 

20.32 

19.89 

> 18.56 

> 

14.35 

> 

14.55 

> 14.15 









Abell 1689 











ID 



9475 

^625 

*775 

•^850 

A 10 

Hl60 


3.6pm 

4.5pm 

5.8pm 

8.0pm 

24pm 

250pm 

350/im 

500/im 

Al 



26.61 

25.78 

24.79 

25.36 

24.29 

23.55 


18.25 

18.45 

18.72 

19.10 

16.94 

> 14.55 

> 14.75 

> 14.35 

A2 



27.35 

26.28 

25.84 

25.59 

> 26.75 

24.36 


> 25.35 

> 24.85 

> 23.95 

> 23.45 

> 20.05 

> 14.55 

> 14.75 

> 14.35 


Note. — The u* magnitudes are MAG_AUT0 values given by our SEXTRACTOR photometry. Similarly, BVR c i' z' magnitudes are MAG_AUT0 values listed in the lFurusawa et al.| j2008f) ’s catalogs. 
The Spitzer photometry values are total magnitudes estimated from the SEIP aperture flux densities with the aperture correction, while the HST photometry values (for A1 and A2) are the 
total magnitudes calculated from MAG_APER values of our SEXTRACTOR measurements. The aperture corrections for the Spitzer and HST photometry are provided by the Spitzer and HST websites, 
respectively (see text). The lower limits correspond to the 3 ex levels. The SI and S2 objects are not observed in u* band and JHK bands. 

t The secure photometry cannot be carried out due to the blending with a bright neighboring objects for the B-500 /im bands. See Figure [l5| 

The positions of the optical-NIR counterpart are far from the ALMA source center by ~ l. // 2. 
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Table 8 

Multiple Images at 2: = 2.60 with the ALMA Counterpart 


ID 

(L07) 

(1) 

ID 

(This work) 
(2) 

R.A. 

(J2000) 

(3) 

Deck 

(J2000) 

cobs 

°1.3mm 

(mJy) 

(4) 

ocorr 

'-'1.3mm 

(mJy) 

(5) 

ocorr 

, 1 . 2 mm 

(mJy) 

(6) 

9475 

7625 *775 

(magAB) 

(7) 

2850 

(8) 

5.1 


13:11:29.064 

-1:20:48.64 

0.18 ± 0.06+ 

0.049 ± 0.021 

0.067 ± 0.029 

26.29 

25.53 

25.37 

25.17 

5.7 

5.2 

Al 

13:11:29.224 

-1:20:44.24 

0.30 ± 0.06 

0.015 ± 0.006 

0.020 ± 0.008 

26.61 

25.78 

25.58 

25.36 

22.8 

5.3 


13:11:34.120 

-1:20:20.96 

< 0.28 

< 0.093 

< 0.13 

26.61 

25.85 

25.85 

25.60 

3.0 


Note. — (1) and (2) ID names given in lLimousin et all J2007D (L07) and this work, respectively. (3) R.A. and Decl. in the optical data. (4) 
Peak flux density SExtractor measurement (Section |3.H with the primary beam correction. If no significant 1.32 mm fluxes are detected at the 
optical source positions, three sigma upper limits are presented. (5) Observed flux density at 1.32 mm on the uv-taper image that are measured 
by 2D Gaussian-fitting routine of imfit in CASA with the primary beam correction. The magnification uncertainty (Section ED is included in 
the error. ( 6 ) Source flux density at 1.2 mm with the flux scaling from 1.32 mm to 1.2 mm (Section[3j. (7) Total magnitudes in the HST /ACS 
bands. ( 8 ) Magnification factor at the peak of the optical flux. 

^ Marginal detection of the 3cr level that is below our source selection limit for the cluster data (> 3.9cr; see Section |3.2| ). 


Table 9 

Very Strongly Lensed Source Candidate 


ID 

R.A 

Deck 

cobs 

°1.3mm 

(mJy) 

ocorr 

‘-’1.3mm 

ocorr 

'-' 1 . 2 mm 

(mJy) 

£475 7*625 2775 

2850 M 


(J2000) 

(J2000) 

(mJy) 

(mag AB ) 


(1) 

(2) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

A3 

13:11:26.897 

-1:20:32.02 

0.24d=0.06 

0.0015 ± 0.0007 

0.0020 ± 0.0010 

> 28.56 > 28.33 > 28.24 

> 27.22 159.5 


Note. — (1) — (7) correspond to (2) — ( 8 ) in Table [ 8 ] These values are for the A3a component (Figure urn . Because A3a and A3b are 
confused in the uv-taper image, these fluxes are obtained on the image without uv-taper. 
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Figure 9. False-color image of Al. The used images and color 
assignment are the same as Figure [I] The center of the cyan circles 
pinpoint the positions of the 5.1 and 5.2 sources, which are two of 
the multiple images o f the galaxy at 2: spe c = 2.60 (Limousi n et al.l 
120071 ; ICoe et aLll2010t) . The radius of the cyan circle is l'/O that 
corresponds to the search radius for the counterpart identifications. 
The red contours show the ALMA 1.3 mm emission at the 2, 3, 4, 
and 5 cr levels. The white curves represent the lines for the lensing 
magnification of \± = 50. 


and i'. The coefficients of x\ and X 2 are 0.314 and 0.207, 
respectively. W ith these p ho tomet ric relations and the 
selection colors dSteidel et al.ll2004r) . the color criteria of 
BX objects are re-written as 

BV —R c i' > -0.2, 

U -BV > BV-R c i' + 0.23, 

BV —R c i' < 0.2(u — BV) + 0.4, 

U -BV < BV-R c i' + 1.0. (12) 

Similarly, the criteria of BM objects are 

BV-R c i > -0.2, 

U -BV > -0.1, 

BV —R c i' < -0.382(17 — BV) + 0.853, 

BV —R c i' < 0.70(17 — BV) + 0.280, 

U -BV > BV-R c i + 0.23. (13) 


We compare these color criteria with our optical- 
NIR counterparts. Among the 17 optical-NIR counter¬ 
parts, 10 are observed with the photometric system of 
u*BVR c i' necessary for the comparison. Removing the 2 
obvious low -2 objects of S4 and S6, we use 8 optical-NIR 



Figure 10. Very strongly lensed source candidate, A3, found 
in A1689. The red contours show the ALMA 1.3 mm emission 
at the 2, 3, 4 cr levels. The white curves denote the critical lines 
for a 2 = 3 source that are estimated with our mass model. Top : 
ALMA continuum image of A3. The cyan tracks show the posi¬ 
tions of multiple images for lensed sources at z = 2, 3, 4, 5, and 
6 that are predicted with our mass model. The ALMA multiple 
image positions (red contours; A3a and A3b) are well reproduced 
by a lensed source at z ~ 3. Bottom : False-color image of HST 
optical data showing the same region as the top panel. The used 
images and color assignment are the same as Figure Q] No optical 
counterparts of A3 are identified. 


counterparts. Figure [TT| displays the two color diagram 
of these 8 sources with the color selection criteria. Six 
sources fall in or near the selection windows of z ~ 2 
BX/BM or z ~ 3 LBGs, where we include Sll source 
at the border of the selection window. The S9 source 
is placed at the upper right corner of Figure [FT] These 
colo rs of the S9 source are similar to z > 4 LBGs (see, 
e.g., ISteidel et al.ll200§l . It is likely that S9 source is a 
LBG at z > 4. On the other hand, S3 clearly escapes 
from the selection windows. We thus find that 7 out of 8 
optical-NIR counterparts are have colors consistent with 
either z ~ 2 BX/BM or z > 3 LBGs. Including the two 
low -z objects of S4 and S6 sources, 70% (= 7/[8 + 2]) of 
the optical-NIR counterparts meet the LBG (BX/BM) 
color criteria. We then compare magnitudes of the 10 
optical-NIR counterparts with the known z ~ 2 BX/BM 
and z > 3 LBGs, omitting u* magnitudes that would 
include a Lyman break whose central wavelength is very 
sensitive to a source redshift. We find that all of the 10 
optical-NIR sources have R c magnitudes in the range of 
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Figure 11. Two-color diagram for our optical-NIR counterparts. 
The cyan and magenta-shaded regions are the BX and BM selection 
windo ws that are def ined in the BV — R c i' and u* — BV color 
plane dLv et al.N201U) . The BX and BM galaxies reside at 2 ~ 
2 — 2.5 and z ~ 1.5 — 2.0, respectively. The yellow shades indicate 
the approximate re gion of the selection window for z ~ 3 LBGs 
dSteidel et al.l2003D , where we assume that the original photometry 
of Uni G and R bands are roughly the same as those of our u *, B , 
and R c bands, respectively. 


Figure 12. Two-color diagram of the BzK galaxy selection in the 
SXDS photometry system The red circles denote the 7 optical-NIR 
counterparts. Star-forming galaxies at 2: ~ 1.4 — 2.5, called sBzK 
galaxies, are selected in the region beyond the solid line defined 
by (z — K ) — (B — z) > —0.2. The pBzK galaxies, passively 

evolving galaxies at 2: ~ 1.4 — 2.5, are selected in the triangular 
region defined with the sold and dashed (z — K > 2.5) lines at 
the upper right corn er. Galactic stars fall in the region b elow the 
dotted line (See also lDaddi et all1 20041 : TYuma et al.ll2012h . 


23 — 27 mag (Tabic [TJ) tha t is typical brightness of z ~ 2 
BX/BM and z > 3 LBGs dSteidel et alJl2(Ml200l . We 
thus conclude that a majority (~ 70%) of the 10 optical- 
NIR counterparts belong to LBG (BX/BM) population. 

6.3.2. Comparison with the BzK Populations 

In contrast with BX/BM galaxies and LBGs, the BzK 
galaxy population includes not only blue star-forming 
galaxies, but also moderately reddened dust-poor star¬ 
forming galaxies (a.k.a s BzK) and pa ssive galaxies (a.k.a 
pBzK) at 1.4 < z < 2.5 dDaddi et al.ll2004j) . We compare 
the color criteria of the BzK population with our optical- 
NIR counterparts. There are 10 out of the 17 optical-NIR 
counterparts that are located in the SXDS region with 
the photometric measurements of 13, z and K bands 
(Table 0E3 Removing the 2 obvious low-z objects of 
S4 and S6, we use 8 optical-NIR counterparts for the 
comparison. 

Because the photometry system of ours (SXDS) is dif - 
ferent from that of the BzK studies dDaddi et alil2004l f. 
we again apply corrections to the color sel ection crite¬ 
ria for BzK galaxies, which are suggested by I Yuma et al.l 

mm, 

(B — z)oaddi+04 = (B — z)sXDS + 0.3, 

(z — -K)Daddi+04 = (z — K) SXDS + 0.1. (14) 

Figure [12] presents the B — z and z — K colors of the 
8 optical-NIR counterparts. We find that 5 out of the 8 

10 Note that the A1689 region misses K band photometry crit¬ 
ical for the BzK selection. 


optical-NIR counterparts meet the color criteria of sBzK, 
while the rest of the three sources, S3, S9 and S10, have 
colors different from sBzK or pBzK galaxies whose red- 
shift range is z ~ 1.4 — 2.5. Because S3 also shows the 
colors different from those of BX/BM galaxies and LBGs 
('Section 16.3.II) . S3 would reside at z < 1.4 that is out of 
the BzK redshift window. Although the colors of S9 and 
S10 are similar to those of stars, in Section 16.3.11 the 
BX/BM galaxy and LBG selections suggest these two 
sources have colors consistent with z > 3 galaxies. Thus, 
S9 and S10 would be galaxies at z > 3 that escape from 
the BzK redshift window. In this way, we confirm that 
the results of the BzK selection are consistent with those 
of BX/BM galaxies and LBGs. In summary, including 
the two low-z objects of S4 and S 6 sources, we find that 
50%(=5/[8+2]) of the optical-NIR counterparts are sBzK 
galaxies at z ~ 1.4—2.5. We then compare magnitudes of 
these optical-NIR counterparts (Table [TJ, and find that 
the optical and NIR magnitudes a re comparable with the 
sBzK galaxies dDaddi et alJl2004) . 

6.3.3. Comparison with the AGN Population 

Six counterparts of SI, S 6 , Sll, S12, S15, and Al are 
detected in the mid IR of 8 and 24 pm bands. Because 
mid-IR bright objects are prominent candidates of AGN 
that have strong hot-dust emission, we compare colors 
of the six mid-IR bright counterparts and a typical AGN 
(Mrk 231) as well as a starburst (Arp 220). Figure [T3l 
shows the colors of the five mid-IR bright counterparts, 
Mrk 231, and Arp 220 in the diagnostic color diagram 
of * 58 . 0 /im/* 54 . 5 /xm vs. ' 8 * 24 //m/‘Ss.o/ym that separates AGN 
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Figure 13. Color-color diagram for the AGN diagnostics. The 
red circles and upper limits are our five mid-IR bright counterparts. 
The dotted and dashed lines present the color tracks of Mrk 231 
(AGN) and Arp 220 (SB), respectively, that are redshifted from 
z ~ 0.5 to 3.0. The numbers labeled to the lines represent the 
redshifts. 


and starburst (SB) populations (llvison et al.ll2004l . In 
Figure [TUI the dashed and dotted lines indicate the color 
tracks of Mrk 231 and Arp 220 in the re dshift range 
of 0.5 — 3.0 that is taken from Figure 4 of llvison et al.l 
(|2004f) . All of our five mid-IR bright counterparts are 
placed near the color track of Arp 220, and far from the 
one of Mrk 231, which suggest that these six counterparts 
are not AGN but starbursts. 


6.3.4. IRX-/3 Relation 

To evaluate the SED properties quantitatively, we 
calculate the infrared-to-UV luminosity ratio IRX(= 
Air/Ai 60 o) and the rest-frame UV continuum slope /3 for 
6 faint ALMA sources of S4, S6, SI 1, S12, S15, and A1 
that have spectroscopic or photometric redshifts. Here 
Lir is the integrated infrared flux density defined in a 
range of rest-frame 8-1000 /im. The /3 values are esti¬ 
mated with two broadband magnitudes, toi and m 2 ' 


P = 


m 1 — TO2 
2.51og(Ai/A?) 


(15) 


where A^ and A^ are the central wavelengths of the two 
broadband filters. We use (mi, m 2 ) = (if*, B), (u*, B ), 
(V, i'), ( V , i'), ( V , i'), (r 625 , zsso) for S4, S6, Sll, 
S12, S15, and Al, respectively. These IRX and (3 values 
are shown in Figure 1141 together with those of nearby 
(U)LIRGs and optically selected star-forming galaxies 
at 2 ~ 2. In Figure Ull we also plot the 5.1 source of 
the multiple image (Table [8]) that is marginally detected 
(3cr). For the lensed sources of Al and 5.1, we assume 
that the positional offsets between the submm and op¬ 
tical sources are real, and use the lensing magnification 
values estimated at the submm (optical) source positions 
to derive the intrinsic submm (optical) luminosities. 


In Figure [bl] all of the six faint ALMA sources fall in 
the region nearly-enclosed by the Calzetti and SMC law 
curves. The six faint ALMA sources are placed clearly 
below the (U)LIRGs that have large IRX values at a 
given /3. The clear difference between the faint ALMA 
sources and (U)LIRGs suggests that the faint ALMA 
sources are not miniature (U)LIRGs with infrared lumi¬ 
nosities fainter than (U)LIRGs by 1 — 2 order(s) of magni¬ 
tude. Note that these results are based on the 59% of the 
faint ALMA sources with the optical-NIR counterparts. 
The rest of the 41% of the faint ALMA sources that 
have no optical-NIR counterparts could have IRX val¬ 
ues as high as those of (U)LIRGs, because these sources 
have low T 1600 values that are suggested by the fact of 
no optical-NIR counterparts. 

These 6 faint ALMA sources are composed of 4 high- 
z sources (Sll, S12, S15, and Al) and 2 low-z sources 
(S4 and S6), where we define high -2 and low -2 by z > 1 
and z < 1, respectively. Four high-z sources have the 
IRX-/3 relation similar to those of optic ally selected star¬ 
forming galaxies at 2 = 2 dReddv et al.ll2012lf . One low -z 
source of S4 is also placed at the region of z = 2 opti¬ 
cally selected galaxies in Figure[]31 while the other low -2 
source of S6 has large IRX and /3 values. It implies that 
majority of high-z sources are similar to z = 2 optically 
selected galaxies, but that low-z sources would be widely 
distributed in the IRX-/3 plot. 

Figure Q3] shows that the multiple images of 5.2 (i.e. 
Al) and 5.1 have similar IRX-/? values. This fact con¬ 
firms that these sources are multiple images of a single 
high-z source. 

6.3.5. What are the Optical-NIR Counterparts of the Faint 
ALMA Sources? 

We summarize the results of Sections I6.3.HIU73.3I in 
Table [TUI We find, in Sections 16.3.11 and 16.3.21 that 7 
out of the 10 optical-NIR counterparts (including the 
low-z galaxies) meet the color criteria of BX/BM LBG 
and sBzK galaxy populations (see Table [TUI) , and that 
the magnitudes of the counterparts are consistent with 
these populations. The rest of the three counterparts 
are low-z objects (2 objects) and probably the fore¬ 
ground/background object of the chance projection (1 
object; Section [63|). These results suggest that not all 
but the majority of the optical-NIR counterparts of our 
faint ALMA sources belong to the sBzKs and BX/BM 
LBGs populations. Moreover, the optical-NIR counter¬ 
parts show neither AGN nor (U)LIRG signatures (Sec¬ 
tions 16.3.31 and 16. 3. 41) . 

We thus conclude that the majority of our faint ALMA 
sources with a detectable optical continuum (~ 25 mag) 
are UV bright star-forming galaxies of BX/BM LBGs 
and sBzKs with no AGN. This conclusion is supported by 
the clustering analysis results in Section [5] that presents 
a weak clustering signal of b g < 3.5 comparable to those 
of LBGs and sBzKs. 

Note that we find the optical-NIR counterparts in 59% 
of our faint ALMA sources fSection lfi.lD . In other words, 
we understand that roughly about a half of the faint 
ALMA sources are LBGs and sBzKs. For the remaining 
half of the faint ALMA sources, the physical origin is 
unclear. Our study identifies a new issue of faint ALMA 
sources that should be addressed in the future studies. 
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Table 10 

Properties of our Optical-NIR Counterparts 


ID 

R.A. 

Decl. 

•^obs 

gobs 

•^obs 

(mjy) 

ocorr 
^ 1.2mm 

•^photo 

^spec 

sBzK 

pBzK 

BX/BM+LBG 

AGN 


(J2000) 

(J2000) 

(mm) 

(mjy) 










(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

SI 

34.612686 

-4.583063 

1.22 

0.36 ± 0.07 

0.43 ± 0.09 






N 

S2 

34.613739 

-4.588505 

1.22 

0.31 ± 0.09 

0.31 ± 0.10 







S3 

34.437805 

-5.494339 

1.03 

0.056 ± 0.015 

0.032 ± 0.010 



N 

N 

N 


S4 

34.446648 

-5.007715 

1.26 

0.047 ± 0.014 

0.044 ± 0.016 

rj 

u ‘' 1 -0.03 






S5 

34.446007 

-5.008159 

1.26 

0.049 ± 0.014 

0.050 ± 0.017 


Y 

N 

Y 


S6 

34.387825 

-5.220457 

1.25 

0.41 ± 0.11 

0.42 ± 0.13 

0 67"*" 0,01 
u ‘° -0.02 





N 

S7 

34.411190 

-4.745051 

1.23 

0.30 ± 0.08 

0.14 ± 0.09+ 


Y' 

N 

Y' 


S8 a 

34.305553 

-5.067293 

1.30 

0.29 ± 0.07 

0.22 ± 0.09 







S9 

34.194702 

-5.056259 

1.25 

0.43 ± 0.12 

0.39 ± 0.15 



N 

N 

Y' 


S10 

34.197121 

-5.059778 

1.25 

0.31 ± 0.08 

0.37 ± 0.10 



N 

N 

Y' 


SlU 

34.442795 

-4.911066 

1.23 

0.53 ± 0.10 

0.56 ± 0.11 

1 54+ u - us 
1 -°^-0.02 


Y 

N 

Y' 

N 

S12 

34.441357 

-4.912045 

1.23 

0.32 ± 0.09 

0.29 ± 0.09 

1,0 -0.02 


Y 

N 

Y 

N 

S13 

34.747295 

-4.858103 

1.25 

0.37 ± 0.10 

0.26 ± 0.11 






S14 

34.274231 

-4.860201 

1.30 

0.32 ± 0.09 

0.43 ± 0.12 







S15 

34.587593 

-5.318766 

1.25 

0.32 ± 0.09 

0.21 ± 0.10 

1 45”*~ 0 ’ 03 
± -^°-0.02 


Y 

N 

Y 

N 

AR 

197.871590 

-1.345783 

1.32 

0.30 ± 0.06 

0.020 ± 0.008++ 

2.60 




N 

A2 

197.869086 

-1.337143 

1.32 

0.31 ± 0.08 

0.077 ± 0.033+t 








Note. — (1): Wavelength in the observed frame. (2): Peak flux density of the SExtractor measurement (Section 1 .1.1 1 with the primary beam 
correction at the observed wavelength. (3): The best-estimate source flux density at 1.2 mm. The source flux is estimated with the 2D Gaussian- 
fitting routine of imfit in CASA (Section I I It with the primary beam correction and the flux scaling to 1.2 mm (Section[3j. For Al and A2 sources, 
the lensing magnification corrections are al so applied. The lensing magnification factors are estimated at t he ALMA flux peak p ositions of Al and 
A2. (4): Photometric redshift estimated bv lWilliams ct all i20Q9i i. (5): Spectroscopic redshift obtained bv ILimousin et all (1200711 - (6), (7), (8), and 
(9): ”Y” indicates the objects that meet the color selection criteria of sBzK, pBzK, BX/BM, and AGN populations, respectively. ”N” represents the 
sources escaping from the color space of the selection criteria. ”Y'” is presented for the sources having colors consistent with the populations. 

+ The photometry of imfit may be biased, due to the systematic noise. The SExtractor peak flux measurement suggests 0.32 =t 0.08 mjy. 
ft The magnification uncertainty (Section 13.ft is also included in the error. 

“Same as the AS2 source in IHatsukade et all 120151 1 
b Same as the AS4 source in IHatsuhade_et al Jl [2015l l 
c Same as the Chen-5 source in IChen et all 120141) 


7. SUMMARY 

In this paper, we identify 133 faint 1.2-mm continuum 
sources with a flux density of 0.02 — 1 mjy in the fifty 
blank-field and one cluster, Abell 1689, maps taken by 
~ 120 deep ALMA Band 6/7 pointings. This is the large 
faint source sample reaching 0.02 mjy, which is drawn 
from our deep ALMA projects and the complete archive 
search for the deep ALMA data that are open for public 
by 2015 June. Using this large sample of the faint 1.2 mm 
sources, we derive the 1.2 mm number counts and discuss 
the faint-source contribution to the EBL. We then inves¬ 
tigate the physical origin of these faint 1.2 mm sources 
by the statistical technique of clustering analysis and the 
individual-source approach of optical-NIR counterparts. 
The major results of this paper are summarized below. 

1. We derive the 1.2-mm number counts in a flux den¬ 
sity range of 0.02 — 1 mjy. Our number counts are 
consistent with the previous ALMA and single-dish 
lensing results, and largely push the flux density 
limit down to 0.02 mjy. We find that the number 
counts are well represented by the Schechter func¬ 
tion at > 0.02 mjy, combining the faint and bright 
number counts from our and the previous studies, 
respectively. 

2. We estimate the total integrated 1.2 mm flux den¬ 
sity with the derived number counts, and obtain 
22.9/^6 Jy deg -2 down to the flux limit of 

This integrated 1.2 mm flux density corresponds 
to 104^25% of the EBL measured by COBE ob¬ 


servations. Because we reach ~ 100% of the EBL- 
resolved fraction at ~ 0.02 mjy within the errors 
that are mainly originated from the COBE EBL 
estimates, the major 1.2 mm EBL contributors are 
sources with > 0.02 mjy. In other words, 1.2 mm 
sources with a flux density below ~ 0.02 mjy make 
a negligibly small contribution to the EBL. These 
results suggest a possibility that there is a flatten¬ 
ing and/or truncation in the number counts at the 
very faint flux density regime of < 0.02 mjy. 

3. We perform the counts-in-cells analysis for cluster¬ 
ing properties of our faint ALMA sources, assuming 
that these sources reside at z ~ 1 — 4. We obtain 
an upper limit of the galaxy bias b g < 3.5 that is 
smaller than those of massive galaxy populations of 
DRGs and SMGs. This upper limit is comparable 
with the galaxy bias values of UV-bright sBzK and 
LBG populations. 

4. We search for optical-NIR counterparts of our faint 
ALMA sources located in two regions of SXDS 
and Abell 1689 that have deep optical-NIR im¬ 
ages and the various multi-wavelength data. We 
identify 17 optical-NIR counterparts, and find that 
59% of our faint ALMA sources have the counter¬ 
parts with the detectable (m < 28 mag) optical- 
NIR continuum. Neither X-ray nor radio counter¬ 
parts are found in any of our faint ALMA sources, 
and the mid-IR diagnostics indicates that there are 
no AGN activities in these sources. We apply the 
























ALMA Census of Faint 1.2 mm Sources 


25 



2-101234 

P 


Figure 14. Infrared-to-UV luminosity ratio IRX as a function 
of the UV slope /3. The red filled circle is the spectroscopically- 
confirmed A1 (5.2) source at z spe c = 2.60. The red open circles 
denote the faint ALMA sources with the photometric redshifts. 
The open black circle is the 5.1 source of the A1 multiple im¬ 
age system that is marginally detected in the ALMA data (Ta¬ 
ble [8]). The blue diamonds (arrow) represent(s) the measurements 
(the upper limit) of the z ~ 2 optically-selec ted galaxies obtained 
by the stacking analysis (Reddy ct al. 2012) . The black squares 
indic ate nearby fUlLIRGs ([Trent ham et al. 1999; Goldad er et al.l 
120021 1. The solid, dashed, and dotted curves denote the IRX-/3 re¬ 
lation of the extinction curves for the Calzetti law (fMeurer et al.l 
1199911 , the updated Calze tti law dTakeuchi et al.l 120121) . and the 
SMC law (Boucliet ct al. 1985). For the Ca lzetti law curve, we 
shift the original relati on of the Calzetti law (Mcurcr ct al. 1999) 
by +0.24 dex, because IMeurer et al.l d 19991) define their infrared 
luminosities at 40 — 120/un that is different from our definition 
(8 - 1000/un). 

color selection criteria of BX/BM LBGs, sBzKs, 
and pBzKs for our optical-NIR counterparts with 
a photometric data set necessary for the color se¬ 
lections. We find that a majority of these optical- 
NIR counterparts have colors meeting the BX/BM 
LBG and sBzK criteria, and that their magnitudes 
are comparable to the BX/BM LBGs and sBzKs. 
Because 59% of our faint ALMA sources have 
the detectable optical-NIR counterparts, we con¬ 
clude that about a half of faint ALMA sources are 
BX/BM LBGs and sBzKs that are high -z galax¬ 
ies with a star-formation rate and dust extinction 
much smaller than SMGs. Moreover, the IRX-/3 re¬ 
lation plot indicates that the optical-NIR counter¬ 
parts of our faint ALMA sources are clearly differ¬ 
ent from (U)LIRGs, and that faint ALMA sources 
are not miniature (U)LIRGs simply scaled down 
with the infrared brightness. 
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Figure 15. The snapshot images of the optical-NIR counterparts. The red contours indicate the ALMA mm band intensity at 2.4, 3.4, 
4.4, and 5.4 a levels. The cyan circle denotes the l'/O search radius for the counterpart identifications. The wavelengths of the images range 
from 0.4 to 500 /mi, which are shown with the labels in the images. The bottom right panels present the ALMA images with the beam 
sizes (yellow ellipses). In each panel, the size of the image is indicated. North is up, and east is to the left. 
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Fig. 1151 (continued) 
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